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Abstract 

The launch of the gamma-ray telescope Fermi Large Area Telescope (Termi-LAT) 
started a pivotal period in indirect detection of dark matter. By outperforming 
expectations, for the first time a robust and stringent test of the paradigm of 
weakly interacting massive particles (WIMPs) is within reach. In this paper, we 
discuss astrophysical targets for WIMP detection and the challenges they present, 
review the analysis tools which have been employed to tackle these challenges, and 
summarize the status of constraints on and the claimed detections in the WIMP 
parameter space. Methods and results will be discussed in comparison to Imaging 
Air Cherenkov Telescopes. We also provide an outlook on short term and longer 
term developments. Keywords: dark matter - galaxies: dwarf - gamma rays: 
galaxies 
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I. INTRODUCTION 

Throughout the course of the past several decades, observational evidence and theoretical 
models have converged to strongly suggest the presence of a dark component of the matter 
content of the Universe, most likely consisting of one or several new particles of Nature. 
Identifying the nature of dark matter is a primary topic of study in modern science. 

A variety of experiments around the world are continually improving their sensitivities, 
ruling out theoretically well-motivated regimes of parameter space for popular extensions 
to the standard model of particle physics that include Supersymmetry and Universal Extra 


Dimensions (Agnese et al, 2014 Akerib et al, 2014; Aprile et al, 2012). As dark matter 


detection experiments continue to improve in sensitivity, new experimental and theoretical 
challenges will continue to be confronted. 

Of particular interest are dark matter particles with weak interactions, labeled as WIMPs 
(Weakly-Interacting Massive Particles). We will here briefly discuss the reason for the 


wimp’s popularity, more in-depth discussions can be found for example in Baer et al 
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(2014) 

Bergstrom 

(2012) 

Bertone, Hooper, and Silk 

(2005) 

Jungman, Kamionkowski, and 

Griest 

(1996 
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At high temperatures (T » mwiMP, where mwiMP is the WIMP mass), particles could 
be thermally created and destroyed, implying a distribution (the Boltzmann law) for 
their state density. As the temperature decreases, the density is exponentially suppressed 
(oc exp —m-wiMp/T) with temperature. Chemical equilibrium is left when the temperature is 
no longer high-enough to pair-create WIMPs, at which point their number density decreases. 
When the WIMP mean free path is comparable to the Hubble distance, the particles also 
cease to annihilate and leave thermal equilibrium, commonly referred to as “freeze-out”. At 
this point, the co-moving density remains constant. Solving for the Boltzmann equation, one 
hnds that the temperature for which the freeze-out occurs is about 5% of the WIMP’s mass, 
thus the density becomes constant (so called “relic density”) when the particles are already 
essentially non-relativistic. Consequently, considering only s-wave annihilation, the relic 
density depends only on the total annihilation cross-section and the velocity distribution: 


Q 


WIMP 




3.1 ■ 10-27 


( 1 ) 


< a a\v\ > 

where the average is taken over velocities and angles. The scale for weak interaction strength 
(~ tt^/’^wiMp) implies that < aA\v\ >~ IQ-^^ cm^s”^, where the WIMP mass is taken to 
be 100 GeV. The resulting relic density for such a particle would be within a factor 3 
of the measured relic density of GwiMph^ = 0.1199 ± 0.0027 provided by the most recent 


PLANCK results (Ade et al ., 2015). This remarkable coincidence is sometimes referred to as 
the “WIMP miracle.” A more careful calculation for the cross-section of thermal WIMPs is 


presented in Steigman, Dasgupta, and Beacom (2012), but does not alter the main argument. 

There are a variety of WIMP search methods ongoing, and for dark matter particles more 
generally. They are broadly classihed as follows: 


• Indirect searches measure the annihilation and/or decay products of dark matter from 
different astronomical environments in the Universe. They specihcally limit the rate at 
which dark matter particles annihilate or decay into standard model particles. Many indi¬ 
rect searches are now underway to detect the remnants of dark matter particles that have 
annihilated into standard model hnal states. The most readily accessible standard model 
particles are photons, in particular high-energy gamma rays. Neutrinos, antiprotons, and 
positrons may also be detected by modern experiments. The annihilation cross section 
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probed by indirect searches is most closely related to the process that sets the dark mat¬ 
ter abundance in the early universe, assuming that the dark matter was once in thermal 
equilibrium. 

• Direct searches measure the scattering of dark matter off of nuclei in low background 
underground detectors, where the collision of dark matter is deduced through energy input 
into particles in the detector. Though the WIMP-nucleus scattering cross section is not as 
simply connected to the dark matter relic abundance as is the annihilation cross section, 
many modern experiments are now probing well-motivated theoretical models, as the generic 
WIMP cross section is determined by its de Broglie wavelength and the weak interaction 
scale. 

• Collider searches measure the production of dark matter through the collision of high 
energy standard model particles, such as protons and electrons. New stable dark matter 
particles produced are identified either by initial state radiation or through the production 
of quarks and gluons, which eventually decay down into the lightest stable particle in the 
spectrum. 

Now that we are in an era in which several different experiments are able to test the¬ 
oretically well-motivated particle dark matter models, it is timely to examine the specihc 
contributions that have been made by the most important experiments. This review article 
will focus exclusively on the recent advances in indirect dark matter searches, and within 
this subject concentrate exclusively on indirect dark matter searches using gamma rays. 
Now that we have had an influx of data over the past few years from these experiments, it 
is important to assimilate this information into the broader context of particle dark matter 
searches. 

Of particular interest are the results that have been obtained by the Fermi Large Area 
Telescope {Fermi-LAT), launched in June 2008 into low Earth orbit. It was designed as 
the successor of the successful EGRET mission, with an order of magnitude better energy 
and angular resolution than EGRET. Fermi consists of two experiments, the Large Area 
Telescope (LAT), which is the primary instrument aboard Fermi, and the gamma-ray burst 
monitor. The LAT is sensitive to photons in the range of approximately 20 MeV — 300 GeV. 
The LAT is an imaging, wide field-of-view pair conversion telescope that measures elec¬ 
tron and position tracks that result from the pair conversion of an incident high-energy 
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gamma-ray in converter foils. The energy resolution over the energy range of interest is 
approximately 10%, and the effective area of the LAT is approximately 10^ cm^. 

In this article we review the constraints on particle dark matter that have been obtained 
from the hrst hve years of the Fermi-LAT and compare them with those obtained by Imag¬ 
ing Air Cherenkov Telescopes (lACT). This article has two main focuses. First, to review 
the experimental methods that are used to analyze Fermi-LAT and lACT. Second, to re¬ 
view the impact that uncertainties in astrophysical dark matter distributions have on the 
interpretation of limits on the dark matter particle mass and cross section. For a recent 
analysis of particle dark matter models that are probed by indirect detection experiments, 
see e.g. 


Bergstrom (2012) 


In our discussion of the experimental methods that are used to analyze the data, we 
review in detail the statistical techniques that have been developed, how they are applied 
to different sources in the universe, and we highlight which of these sources provide the best 
detection prospects. 

In our discussion of astrophysical systematics, we review methods to determine dark 
matter distributions in all of the various sources that have been studied by FermALAT. We 
discuss not only how the mass distribution of the Milky Way, dwarf galaxies, and galaxy 
clusters are determined, but also look forward to how future astronomical surveys and 
measurements will improve the understanding of dark matter distributions. We discuss 
the connection that will continue to be strengthened between astroparticle experiments and 
larger scale astronomical surveys. 

This review is organized as follows. In Section |TT] we discuss the different astrophysical 
targets that are used for indirect detection, and the pros and cons of each. In Section III, we 


review different types of gamma-ray instruments, and discuss the analysis challenges that 
they face. In Section [IV| we review the current status of indirect dark matter searches with 
gamma rays. In Section |V], we discuss what future gamma-ray experiments will bring to the 


held of indirect dark matter detection, and in Section VI we present conclusions. 


II. TARGETS: PROMISES AND CHALLENGES 


Astrophysical uncertainties have long presented a systematic uncertainty in searches for 
new physics. Focusing on gamma-ray searches for dark matter, there are two broad types of 
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systematic uncertainties that must be accounted for: uncertainties in modeling of dark mat¬ 
ter distributions, and uncertainties in gamma-ray emission from non-dark matter sources. 
In this section, we primarily discuss the contribution of the former type of uncertainty, 
and defer to discussion of the latter type of uncertainty when we present the results of the 
gamma-ray observations in the following sections. 


Within the context of the A-Cold Dark Matter (ACDM) paradigm (Frenk and White 


2012), observed galaxies are formed within dark matter halos. These halos are formed 


through a sequence of mergers of lower mass halos and a process of smooth accretion onto 
the halo. As a result, dark matter halos are complex systems that are not totally “smooth,” 
but rather have features in their phase space in the form of subhalos, or substructure, and 


tidal debris that reflects their interaction within the larger halo (Diemand et al., 2008 


Vogelsberger et al, 2009). Dark matter halos are not predicted to be perfectly spherical. 


but rather retain shapes that reflect their formation and interaction histories (Allgood et al 


2006). 


Though not spherical, it is often convenient to express the density distribution of dark 
matter halos in a spherically averaged form, for example as 

Ps 


p{r) = 


( 2 ) 


(r/rs)'i'(l -h (r/r5)^)h-DA' 

Here ps and are scale density and scale radius parameters, respectively, h represents the 
turnover from the asymptotic power-law slope in the inner regime, 7 , and the asymptotic 
power-law slope in the outer regime, c. The set of parameters ( 7 , 6 , 0 ) = (1,1,3) give the 
well-known Navarro-Frenk-White prohle (NFW, Navarro, Frenk, and White, [1997 ). 

Though an NFW prohle is now a robust prediction of dark matter-only simulations of 
halos, as we discuss below in this section there is debate as to how well this prohle describes 
observations of the dark matter halos over a wide range of mass scales. In spite of this 
debate. Equation provides a useful starting point for phenomenological studies of dark 
matter halos, and will be referred to throughout the course of this section and review. 

With a well-motivated model for dark matter halos in place, we can make predictions for 
gamma-ray huxes, and the corresponding uncertainties, from diherent astrophysical targets. 
We in particular focus on the uncertainties in the dark matter properties of each of these 


systems, and thus the observational uncertainties on the quantity. 


fAO 


J(AD) = - 


r-A 


sin TdT 


p^ 


(3) 
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where ^ represents the angular separation from the center of the halo, which is typically 
deduced from the position of the center of the observed galaxy. Here D is the distance to 
the center of the galaxy, so that — 2iD cos d'. The upper and lower boundaries 

to the integral are i± = D cos T ± y/r^ — D'^ sin^ T, where r* is the tidal radius of the dark 
matter halo. Equation]^ will be referred to here as the “J-factor.” This sets the dependence 
of the annihilation signal on the dark matter distribution in any astrophysical system. As 
we highlight in this section, the determination of the density prohle from observation is 
not straightforward, and thus the uncertainty in the J-factor translates into a systematic 
uncertainty in the gamma-ray flux determination. 

We begin by examining the dark matter distribution in the Milky Way galaxy, includ¬ 
ing the Galactic center. We then move on to discuss determination of the dark matter 
distributions in the dwarf spheroidals of the Milky Way, which are probably the simplest 
astronomical systems from which the dark matter distributions can be derived. We then 
follow up with a discussion on dark matter substructures. Galaxy clusters, the Galactic and 
the extragalactic dark matter distributions. 


A. Galactic Center 


Most likely, the largest flux of gamma rays from dark matter annihilation comes from 
the Galactic center because of its high concentration of dark matter and close proximity. 
However, a drawback of the Galactic center as a target is that there is a substantial pop¬ 
ulation of gamma-ray sources and there is diffuse emission from cosmic rays that must be 
well-understood in order to extract the dark matter signal. An additional drawback of the 
Galactic center as a target for gamma-ray and dark matter studies is that it has proven 
challenging to extract the dark matter distribution in the Galactic center. In the central 


few parsecs, the stellar mass from the bulge dominates the dynamics (Schodel, Merritt, and 


Eckart, 2009). Even including a weak disk component, the shape of the dark matter distri¬ 


bution is unconstrained, so that it is not possible to tell if the dark matter prohle rises to 


a central cusp-like structure or has a constant density within the bulge region Wocco et al 


2011). A conservative upper bound on the contribution of the dark matter may be set by 


the upper limit deduced from the bulge contribution to the potential. 

From a theoretical perspective, the dark matter distribution near the black hole at the 











center of the Milky Way is modified by the black hole itself (Gondolo and Silk, 1999). 
Adiabatic growth of the central black hole may steepen the central cusp of dark matter. 
Assuming an initial dark matter prohle in the Galactic center of the form p oc r~"‘^ from 
conservation of mass and angular momentum, the hnal mass prohle scales as p oc where 


A = (9 — 27)/(4 — 7 ) (Ullio, Zhao, and Kamionkowski, 2001). Gondolo and Silk (1999) derive 
a lower bound of ~ 0.24 GeV cm“^ of dark matter near the Galactic center. This result 
relies on assumptions that are not well understood, for example the adiabatic nature for 
the growth of the black hole itself, and scattering of dark matter off of stars in the central 


nuclear star cluster (Bertone and Merritt, 2005). It is also possible for dark matter to 
interact with the central nuclear star cluster and form a prohle that is shallower than the 


NFW prohle (Gnedin and Primack, 2004). 

Numerical simulations have examined the ehect of baryonic physics on the structure 
of Milky Way-mass dark matter halos. Though the inclusion of baryons leads to many 
challenges, numerical simulations have continued to make progress. Simulations of single 
halos, in combination with analytic models, show that the central densities of dark matter 
halos become less steep than those found in pure N-body simulations because the baryons 


induce repeated epochs of feedback due to star formation activity (Pontzen and Governato 


2012 ). 


The lack of consensus from both the observational and theoretical side on the nature of 
the Milky Way dark matter density prohle provides a signihcant systematic uncertainty that 
must be accounted for in gamma-ray searches for dark matter. In fact it is likely that, from 


a pure astronomical perspective, this lack of consensus will remain for some time. locco. 


Pato, and Bertone (2015) have recently compiled known data to analyze the dark matter 


content within the Solar circle, but are insensitive to the dark matter density prohle. Their 
claim that this estimate constitutes the hrst evidence for dark matter has however triggered 


some debate(McGaugh et al. 2015). 

Kinematic data from the GAIA satellite is expected to somewhat improve on the system¬ 


atic uncertainty in the measurement of the local dark matter density (Bovy and Rix, 2013), 
though it will be more difficult to use this data to determine the shape of the prohle near 
the Galactic center. 
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B. Milky Way mass profile 


The discussion above focused on the central density profile of the Milky Way, i.e. from 
the Solar circle towards the Galactic center. Going in the other direction, the density profile 
of the Milky Way as measured from the Solar radius out to the virial radius is important 
for determinations of diffuse gamma-ray emission. However, because of our position within 
the Milky Way’s dark matter halo, it is more difficult to determine its dark matter density 
prohle than it is for many external galaxies. 

The best measurements of the integrated mass of the Milky Way come from spectroscopy 
of stars in the outer region of the dark matter halo. The Sloan Digital Sky Survey (SDSS) 
survey has measured the dark matter mass prohle of the Milky Way using the kinematics 


of a large sample of Blue Horizontal Branch (BHB) stars (Xue et al ., 2008). BHB stars are 
important tracers of the Milky Way mass because they are both intrinsically bright and have 
accurate distance measurements. Due to the average radius of the BHB stars, they provide 
the best constraint on the Milky Way mass within a radius of about 60 kpc, measuring a 


mass of ~ 4 X IO^^Mq. This is consistent with the results of independent analyses (Deason 
et al ., 2012), and with estimates at larger radii, which hnd a mass of ~ 7 x IO^^Mq enclosed 


within 80 kpc ( Gnedin et al.\ 2010). Several authors have used samples of bright BHB 
stars to extrapolate and determine the total mass of the Milky Way. They have reported a 
total mass of the Milky Way ranging anywhere from 0.5 — 2.5 x IO^^Mq (|Battagha et al. 


2005 Busha et al., 2011; Xue et al., 2008), with the variation depending on the exact 


analysis method and the sample of stars used. These results can be compared to updated 
implementation of the timing argument, which implies a total Milky Way mass of ~ 2 x 10^^ 


Mq (Li and White, 2008). 

Distant satellite galaxies can also be used as tracers of the mass distribution of the Milky 


Way (Little and Tremaine, 1987). For these measurements, there are two systematic un¬ 
certainty that are particularly important. The first is the due to the uncertain density 
distribution of the satellites in the Milky Way; this is mainly due to the small number of 
satellites that are currently known, about a couple dozen. The second, and perhaps proba¬ 
bly the most important systematic uncertainty, involves understanding whether the distant 
satellite Leo I is bound to the Galaxy. Because Leo I is moving at a high Galactocentric 
velocity and is at a Galactocentric radius of 260 kpc, it is not yet clear whether this is an 
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outlier that is bound to the halo or if it is unbound and is on its hrst pass through the 
Galaxy. Several recent proper motions do in fact seem to indicate that it may be bound to 


the Galaxy (Boylan-Kolchin et al ., 2013). The Milky Way mass measurements from both the 


BHB stars and the satellite galaxies are in good agreement with the recent measurements of 


the Milky Way escape velocity using a local sample of high velocity stars (Piffl et al., 2014) 


and constraints using the Sagittarius tidal stream (Gibbons, Belokurov, and Evans, 2014). 


Similar to the case of the Galactic center, the factor of a few uncertainty that still lingers 
in the measurement of the mass of the Milky Way’s dark matter halo provides a systematic 
uncertainty in the predicted gamma-ray emission from dark matter annihilation. 


C. Dwarf Spheroidal galaxies 

There are now approximately two dozen satellite galaxies of the Milky Way that are 
classihed as dwarf Spheroidal (dSphs). The dSphs are conveniently classihed according to the 
period of discovery. The hrst nine dSphs that were discovered before the turn of the century 
have come to be classihed as “classical” satellites, while those discovered in the era of the 
SDSS are largely referred to as “ultra-faint” satellites. The dSphs range in Galactocentric 
distance from approximately 15 — 250 kpc, and their overall distance distribution in the 
Galactic halo is much more extended than the more centrally-concentrated globular cluster 
population. Nearly all of the dSphs are devoid of gas up to the present observational limits. 


For a more thorough discussion of astrophysical aspects of dSphs see McGonnachie (2012a); 


Strigari 

(2013) 

Walker 

(2013) 


Dwarf spheroidals provide excellent targets for gamma-ray searches for WIMPs for several 
reasons. First, theoretically there is expected to be no gamma-ray point sources and no 
intrinsic dihuse emission associated with them. Second, their dark matter distributions are 
directly derived from the stellar kinematics. Third, the boost factor from dark matter sub¬ 
structure is predicted to be negligible in these systems, so the interpretation of the limits 
(or detections) from them is much more straightforward than it is for other astrophysical 
targets (see more detailed discussion below on boost from halo substructure). 

Focusing specihcally on dSphs, it is convenient to start by assuming that their stellar dis¬ 
tributions and dark matter distributions are spherically-symmetric. It has become standard 
to model the measured line-of-sight velocities of stars in dSphs with the spherical Jeans 


11 





















equation, and thus extract a measure of the dark matter content. Though in the Jeans 
equation the mass prohle is the sum of the contribution from the dark matter and the stars, 
in nearly all cases of interest for gamma-ray studies the stars contribute little to the total 
gravitational potential. The measured line-of-sight velocity of a star in a dSph is a mixture 
of the radial and tangential velocity components, and therefore depends on the intrinsic 
velocity anisotropy of the stars. The anisotropy is traditionally dehned as the following 
parameter 


(Jt 






( 4 ) 


where is the radial velocity dispersion and cr^ is the tangential velocity dispersion. At a 
projected position i?, the line-of-sight velocity dispersion is 

2 f 


(^losiR) = 


is{R) 


fR 


td2 


Ps(Tlr 


. Vr 




zdr. 


( 5 ) 


The three-dimensional radial velocity dispersion is determined from the Jeans equation, 
and the projected stellar density prohle, Is{R), is determined from hts to the stellar number 
counts as a function of projected radius. It is manifestly clear from Equation that there 
is a degeneracy between the mass prohle and the velocity anisotropy of the stars. 

The half-light radii for the classical dSphs are typically of the order few hundreds of 
parsecs, and their luminosities spread over a range of nearly two orders of magnitude, ap¬ 
proximately 10^ —10^ Lq. The photometric prohles as derived from star counts are typically 
consistent with a cored model in projection, followed by a turnover into an exponential fall- 
oh in the outer region of the galaxy where the stellar density blends into the background 
star counts. For the classical dSphs there are hundreds, and in some cases thousands, of 


bright giant stars that have measured velocities to a precision of a few km/s or less (Walker, 


Mateo, and Olszewski, 2009). 


Constraints on the dark matter distribution of dSphs have been specihcally calculated via 
the following procedure. The kinematic data from a dSph, below dehned as Vk, comprises 
of n stars each with a measured line-of-sight velocity n* and uncertainty at a projected 
radius i?*. Dehne a as a set of theoretical parameters that are to be extracted from the data. 
The probability for the data, given the model parameters, is 


VkiniPkla) oc JJ 
1=1 


- 0 . 5 - 






+ 


{R^ 


exp 




( 6 ) 
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where v is the mean velocity. The projected line-of-sight velocity dispersion, aios{R), which 
is calculated from the spherical Jeans equation, depends on the parameters of the dark 
matter density prohle. A standard assumption for the dark matter density prohle is the 
generalized double power law model in Eq. 

With aforementioned improvements in observational data and theoretical modeling of 
dSphs around the Milky Way, it has become timely to search for the presence of NFW-like 
dark matter density prohles predicted by A-CDM. However, in spite of the high-quality 
modern data sets, there has been signihcant debate as to whether the data are unable 
to uniquely specify a model for the dark matter potential. Indeed, the kinematic data 
of several dSphs, when modeled as a single stellar population with spherical Jeans-based 


models, is consistent with both cusped dark matter prohles (Strigari et al ., 2007) and cored 


prohles (Gilmore et al.. 2007; Richardson and Fairbairn, 2013). 

The assumption of the spherical Jeans equation is of course an approximation to the true 
dynamical state of the dSph. Indeed, all of the dSphs are observed to be elongated, with an 


approximate 30% diherence between the length of the major and minor axis (McConnachie 


2012b). Further, the solutions to the spherical Jeans equations need not admit fully self- 


consistent dynamical solutions with a positive dehnite stellar distribution function. Recent 
work has relaxed the assumption of spherical symmetry in the Jeans equations, in most cases 


hnding that the central mass distributions are consistent with the spherical case (Hayashi 


and Chiba, 2012). There has also been recent work dedicated to establishing self-consistent 


distribution function models of the dSphs. For example, it is possible to hnd self-consistent 
solutions in which the orbits of the stars are isotropic, and the stars trace the dark matter 


prohle in the central region (Strigari, Frenk, and White, 2010). Orbit-based models of 
the dSphs are now being developed; these types of models are ideal in the sense that the 


distribution function is obtained in a model-independent manner (Breddels et al., 2013 


Jardel et al ., 2013). All of this modeling indicates that, when considering the observed stars 


as a single stellar population, the mass estimates are in good agreement with those obtained 
through the Jeans equation, and further it is not possible to conclusively establish whether 
there is a dark matter core or cusp within these systems. 

Some dSphs exhibit evidence for more than a single stellar population, i.e. populations 
of stars that formed at different epochs. Two well-known examples are Sculptor and For¬ 
nax, with Sculptor, which is located at an approximate Galactocentric distance of 80 kpc. 
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receiving a significant amount of attention. Phase space modeling of the kinematics of its 
stars, when taken as a single stellar population, indicate that it is possible to find NFW- 


based dark matter prohles with isotropic stellar velocity dispersions (Breddels et al, 2013 


Strigari, Frenk, and White, 2010). Due to high quality information on both the kinematics 


and the metallicity of stars in Sculptor, it is possible to break up the stars into distinct 
populations. In particular, Battaglia et al. ( 2008[ ) (BOS) have shown that there are two dis¬ 
tinct stellar populations; a metal rich population that is centrally-concentrated, and a more 
extended metal-poor population. Using Jeans-based modeling, BOS showed that in order 
for both populations to be embedded into a single dark matter halo, both the metal rich 
and metal poor populations must transition from isotropic in the center to predominantly 
radially-biased orbits in the outer regions. Following up on the BOS analysis, Agnello and 


Evans (2012) apply the projected virial theorem to the Sculptor data, and hnd that it is 


not possible to self-consistently embed both populations into a halo with an NFW density 


prohle. Further, Amorisco and Evans (2012) study Michie-King models for the stellar dis¬ 
tribution function, and find that even though NFW models provide an acceptable fif fo 
the data, in general cored models are preferred. 

The analysis of BOS was then soon followed up upon with an independent and larger 


data set by Walker and Penarrubia (2011) (WPll). With their data set, WPll also present 
evidence for two populations, with seemingly similar velocity dispersion and half-light radii to 
what was obtained by BOS. Though there is no observed evidence for two stellar populations 
in the WPll data, they extract the dispersion, half-light radii, and metallicity of their 
populations using a statistical algorithm. WPll then apply the mass estimator presented 
in Walker et al. (2009) and Wolf et al. (2010) to the two populations, finding that an NFW 


model is ruled out at the 99% c.l. While this evidence for a dark matter core is intriguing, 
and has garnered a substantial amount of attention in the dark matter community, it has 
not been able to stand up to more rigorous modeling of the two populations in Sculptor. 
In particular, it is possible to find a self-consistent stellar distribution function model with 
an NFW dark matter profile that is able to statistically-describe the two populations in 


Sculptor (Strigari, Frenk, and White, 2014). 

The above discussion underscores the difficulty in determining whether dark matter cores 
or cusps exist in dSphs. While the kinematic data is unable to determine the slope of the 
central density of the dark matter in dSphs, it is much more effective at determining the 
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integrated mass within the half-light radius, approximately a few hundred parsecs (Walker 


et al ., 2009 Wolf et al. 


This is weakly dependent on whether there is a central 
core or cusp in the dSph, so long as the log-slope is —(i(logp)/(ilogr < 1.5. This implies 
that the constraints on the mass prohle directly translate into constraints on the J-factor in 
Equation within the same region (Strigari et a/., 2007 Strigari et al. ^ 2008). For a dSph 
at a distance of ~ 50 — 100 kpc, the half-light radius corresponds to less than approximately 
one degree, which is about the angular resolution of the Fermz-LAT over a large energy 
range of interest. This is the region within which the integrated density and the integrated 
density-squared are the best constrained from the kinematic data sets. This implies that 
the assumption of a core or a cusp for the density prohle does not signihcantly affect the 
gamma-ray hux predictions for the Fermi-LAT. 

Several authors have now published an analysis of the dark matter distributions in dSphs 
using the spherical Jeans method and examined their implications for gamma-ray exper¬ 


iments (Charbonnier et al.. 2011; Essig, Sehgal, and Strigari, 2009, Essig et al. 2010 


Geringer-Sameth, Koushiappas, and Walker 2015 Martinez] 20131 Strigari et al ., 20071 Stri- 


gari et al ., 2008). A typical approach is to assume a model for the dark matter prohle, such 


as a generalized double power law model, combined with the likelihood function in Eq. 
Within a Bayesian framework, the model parameters Ps,'f's, a, b, c are then marginalized over 
assuming priors on these parameters. In the literature there have been several approaches 


to handle these priors. Strigari et al. utilized priors from CDM simulations (Strigari et al. 


2007 Strigari et al.. 2008), which ehectively weighted Eq. [^with a function describing the 


relation between ps, Vg derived from CDM simulations. Several authors have considered “un¬ 


informative” priors, equivalent to hat priors on logr^ and logp^ (Charbonnier et al.. 2011 


Essig, Sehgal, and Strigari, 2009 Essig et al. 2010). Martinez introduced a hierarchical 


modeling method that uses a relationship between that mass at the half-light radius and 
the luminosity of a dSph (Martinez, 2013). As shown in Figure for dSphs with well- 
measured kinematics, the J-factors that are derived for each prior are typically consistent 
with one another, though there is a larger spread for dSphs with smaller samples of stars. 
As noted above these results are weakly dependent on whether a cored or cusped central 
density prohle is assumed for the dark matter. 

Figure [^clearly indicates which dSphs are the most interesting targets for indirect dark 
matter detection experiments. The two dSphs with the largest J-factors, Segue 1 and Ursa 


15 
























































































Major II, are ultra-faint satellites with sparse samples of stars associated to them (about 
60 and 20 stars, respectively). Though these dSphs have the largest mean flux, they also 
have the greatest uncertainty due to the small stellar samples. After Segue 1 and Ursa 
Major II, the dSphs with the next largest J-factors are Ursa Minor and Draco, at 66 and 
80 kpc, respectively. These J-factors are determined from samples of hundreds of stars so 
their corresponding uncertainties are much lower than the uncertainties on the J-factors for 
Segue 1 and Ursa Major II. 



FIG. 1 J-factors within 0.5° for dSphs, for three different assumptions for theoretical priors on 


the dark matter halo parameters. Blue points utilize uniform priors in logr^ and logp^ (Abdo 


et al., 2010c), red points uses as priors the relation between Ps,rs that is determined in CDM 


simulations (Ackermann et al., 2011), and black points uses the hierarchical modeling method that 


is introduced in Martinez (Martinez, 2013). Though these calculations assume an NFW profile, 
for J-factors determined within an angular region of 0.5°, the results are weakly-dependent on the 
assumption for the dark matter density profile. 


In sum, because of the substantial theoretical and observational work that has gone into 
understanding the kinematics of dSphs and their underlying dark matter distributions, they 
are a unique target in gamma-ray searches for dark matter. Any possible detection of a 
signal in other sources must be corroborated by a detection in dSphs. As discussed in the 
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sections below, dSphs are targets for several gamma-ray experiments which hope to probe 
the WIMP mass range from a few GeV up to tens of TeV over the next decade. 


D. Substructures 


N-body simulations of Galactic halos predict that approximately 10 — 50% of the dark 
matter mass of the Milky Way is bound up in the form of substructure, or subhalos. Accord¬ 
ing to the AGDM model, some fraction of these subhalos should host the observed dSphs that 
were discussed above. However, there is not a consensus understanding as to what the mass 
and radial distribution is for the subhalos that “light up” with stars- this issue is strongly 


intertwined with the classical missing satellites problem (Kauffmann, White, and Guider- 


doni, 1993 Klypin et al., 1999 Moore et ai, 1999) and the more recent too-big-to-fail issue 


of AGDM (Boylan-Kolchin, Bullock, and Kaplinghat, 2011). Given that subhalos without 


any associated stars have yet to be conclusively detected around the Milky Way, a gamma- 
ray signal from these objects, while intriguing, is still subject to a substantial amount of 
theoretical assumptions. Nonetheless, it is informative to examine what the modern theory 
predicts for this population of subhalos around the Milky Way, and beyond. 

The numerical simulations of Milky Way-mass galaxies are now complete in their measure¬ 
ment of subhalos down to a mass of approximately 10® Mq, corresponding to approximately 
10“® times the total mass of the host halo ( |Diemand et al. 2008; Springel et ai, 2008a). 
The mass function of subhalos may in fact extend down to Earth masses or even below. 


which the simulations are not sensitive to at present (Green, Hofmann, and Schwarz, 2004 


Loeb and Zaldarriaga, 2005 

Martinez et al. 

[2009 

Profumo, Sigurdson, and Kamionkowski 


2006). For for subhalos with mass greater than ~ 10® Mq, the mass function of subhalos 


is a power law that scales as dN/dM oc M“", with a = 1.9. Because 1 < a < 2, subhalos 
at the low mass end of the mass function dominate the distribution by number, while the 
subhalos at the high mass end of the mass function dominate the total mass in substructure. 

The internal density profiles of the dark subhalos is important, and has been a subject 
of theoretical debate. More recent simulations hnd that dark matter halos and subhalos are 


better described by an appropriate shallower Einasto density prohle (Springel et al ., 2008b). 


There are suggestions from numerical simulations that the density prohles of the smallest, ~ 
10“® Mq halos, are steeper than the NEW form (Anderhalden and Diemand, 2013; Ishiyama, 



























































































Makino, and Ebisuzaki, 2010). The gamma-ray signal from WIMP annihilation from small 


dark matter halos is also sensitive to the concentration of low mass halos, where here the 
concentration is dehned as the ratio of the scale radius of the halo to its virial radius. There 
is a substantial amount of uncertainty on the extrapolation of the concentration versus 
halo mass relation down to low halo mass scales; numerical simulations are at present only 


able to measure this relationship down to halo mass of about 10® Mq (Gao et al, 2012 


Springel et al, 2008b). There is at least an order of magnitude uncertainty in the gamma- 


ray signal from dark matter annihilation in small halos due to the unknown extrapolation of 
the concentration versus mass relation to the smallest halos (Martinez et al. [2009; Sanchez- 


Conde and Prada, 2014). 


Dark matter substructures are very much “high risk, high reward” targets for gamma- 
ray searches. While a conhrmed detection of an object that shines only in gamma rays 
would be a spectacular conhrmation of both the ACDM and the WIMP paradigm, for a 
canonical thermal relic scale WIMP cross section predictions for signal detectability vary 
by several orders of magnitude (]Koushiappas, Zentner, and WalkCT 2004 Pieri, Bertone, 


and Branchini, 2008 Springel et al., 2008b). Realistic a priori predictions for a signal 


from dark matter substructure will only likely be available once simulations are able to 
resolve substructure near the Earth mass scale, or it becomes observationally feasible to 
detect substructure at this mass scale via gravitational lensing. Since there are substantial 
hurdles in reaching these scales in both the simulations and observations, probably the best 
these types of searches can hope for at the stage is to hnd a signal in gamma rays from a 
substructure candidate that can be followed up on and found not to be observed at other 
wavelengths. 


E. Galaxy Clusters 

The study of dark matter in galaxy clusters has a long history, dating back to the original 
discovery of Zwicky. The systematic uncertainty in the determination of dark matter mass 
prohles and J-factors in clusters is similar in many ways to the case of dSphs discussed 
above. The dominant uncertainty in the J-factors from clusters arises from two orthogonal 
aspects of astrophysics. First, there is an uncertainty in the empirically-measured cluster 
mass prohles, which are derived from a combination of x-ray temperature prohles and gas 
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kinematics. Second, there is significant uncertainty in the predicted gamma-ray luminosity 
that arises from the contribution of dark matter substructure in the clusters. 

Cluster masses and density profiles are measured through x-ray emission, galaxy dynam¬ 


ics, or gravitational lensing (Allen, Evrard, and Mantz, 2011). For nearby clusters, masses 


are generally derived under the assumption of hydrostatic equilibrium, 

kT{r)r 


M{r) = 


fimpG 


din p{r) ^ dInT(r) 


dlnr 


dlnr 


(7) 


where T{r) is the temperature profile, and p is the mean mass per particle in units of the 
proton mass. The assumption of hydrostatic equilibrium is subject to systematic uncertainty 


depending on the physical state of the cluster. For instance, Nagai, Vikhlinin, and Kravtsov 


(2007) used numerical simulations to show that systematic uncertainties due to non-thermal 
pressure introduce a ~ 10 % systematic uncertainty in the mass determination. 

Furthermore, there are only of order tens of clusters that have measured temperature and 
density profiles allowing the use of Equation and only a handful of nearby clusters have 
dark matter mass profiles constrained by redundant estimates. For the majority of known 
clusters, we are instead left with indirect observational proxies that are calibrated by low- 
redshift clusters. Such standard observational proxies are the average temperature, the mass 


of the hot gas, and the product of these two, the total thermal energy (Vikhlinin et al ., 2009). 
Recent studies have combined all of the aforementioned mass measurement techniques to 


obtain an estimate of the slope of the dark matter density profile in clusters (Newman et al. 


2013). In these studies, the measured mean slope is ( 7 ) ~ 0.50, which is less shallow than 


the central slope of the standard NFW profile. 

Of particular importance for gamma-ray analyses are the clusters that are the appropriate 
combination of the most nearby and the most massive. From their measured mass distri¬ 
butions and known distances, several recent studies have come to the agreement that the 
Fornax, Coma, and Virgo clusters are the brightest source of gamma rays from dark matter 
annihilation (Ando and Nagai, p012; Han et al ., [2012; Pinzke, Pfrommer, and Bergstrom 


2011 ). 


Dark matter-only simulations of clusters of galaxies are only able to resolve dark matter 
substructure with mass 10 “^ times the mass of the host cluster; this is perhaps 10 orders of 
magnitude larger than the minimum mass dark matter subhalo that is predicted in ACDM 
theory. This uncertainty due to the dark matter substructure, along with the fact that only a 
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relatively small number of clusters has been simulated, introduces a significant uncertainty 
in the contribution of cluster substructure to the gamma-ray luminosity that arises from 
dark matter annihilation. 

Nevertheless, a signihcant increase in the gamma-ray signal from clusters is expected 
from the presence of such dark matter substructure. Numerical simulations provide the 
most reliable method for determining the gamma-ray emission from subhalos in clusters. 
From a sample of nine cluster-mass dark matter halos, which have individual particle masses 
of ~ 10® Mq and identify subhalos down to a mass scale of ~ 10^ Mq, Gao et al. (2012) 
directly determine the overall boost factor and the surface density prohle of the substructure 
component for clusters as a function of the cluster mass. Though it relies on an extrapolation, 
the substructure likely implies a substantial increase in the gamma-ray luminosity over the 
smooth component in clusters, perhaps up to several orders of magnitude. 

In sum, clusters represent a unique target for dark matter searches. Like the dSphs, they 
can be localized in space, and their dark matter distributions can be robustly measured 
from astronomical data sets. However, they are different from dSph targets because there is 
expected to be a large contribution from dark matter substructure within them. While this 
is expected to increase the gamma-ray signal, in the case of a null detection it does provide a 
substantial systematic uncertainty when attempting to set upper limits on the annihilation 
cross section. 

Clusters are also distinct from dSphs because there is expected to be a signihcant hux 
of gamma rays due to cosmic ray processes. The gamma-ray luminosity due to cosmic rays 
in clusters is expected to trace the gas density, so it is more centrally-concentrated than 
an expected dark matter annihilation signal, since the dark matter signal is more extended 
because of the emission from subhalos. As we discuss in the sections below, though very 
plausible theoretical predictions show that gamma rays from both cosmic rays and dark 
matter could have been detected by gamma-ray observatories, there has been no conclusive 
signal reported to date. 

F. Cosmological mass function of dark matter halos 

The discussion above has focused on measurements of the dark matter distribution within 
the Milky Way, and nearby identihed dwarf galaxies and galaxy clusters. Gamma-ray 


Gao et al. (2012 
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searches for dark matter are also sensitive to the accumulated emission from dark mat¬ 
ter annihilation in all the dark matter halos that have formed in the universe. Though this 
emission is “unresolved,” it is possible to model given a mass function of dark matter halos. 
It may be possible to deduce this mass function from the observed luminosity function of 
galaxies, though this method is hindered by uncertainties in the mapping of galaxy lumi- 


nosity to dark matter halo mass ( 

Behroozi, Gonroy, and Wechsler, 

2010 

Guo et al. 

2010 

Kravtsov, Vikhlinin, and Meshscheryakov 

2014; 

Moster et al. 

2010 

)• 


A more robust estimate of the dark matter halo mass function comes from large scale 
cosmological simulations. The two largest volume cosmological simulations to date provide 
a statistically-complete sample of dark matter halos down to a maximum circular velocity 


of about 50 km/s, or a mass of approximately 10® Mq (Boylan-Kolchin et al. , 2009 Klypin, 


Trujillo-Gomez, and Primack, 2011 Springel, 2005). These studies hnd that the “halo 


multiplicity function,” which is derived from the halo mass function, appears to have a 
near-universal form at all redshift. This implies that it is possible to compute the halo 
multiplicity function at any redshift from well-measured cosmological parameters. 

As discussed above in the context of galaxy clusters, dark matter annihilation signals 
are sensitive to substructure within dark matter halos. This implies that it is important to 
determine not only the mass function of dark matter halos, but also the mass function of 
dark matter substructure. Both the shape and the scatter in the subhalo mass function in 
simulations is well-characterized down to subhalo masses at least three orders of magnitude 
less than the his halo mass. Semi-analytical models are now able to reproduce the trends 


observed in the simulations (van den Bosch and Jiang, 2014). 

Given the uncertainties associated with the distribution of dark matter substructure in 
halos, in particular in the extrapolation of substructure down to the smallest mass scales 
below the resolution limit of simulations, theoretical predictions for the diffuse gamma-ray 


background from dark matter halos is difficult to precisely pinpoint (Sefusatti et al., 2014 


Ullio et al., 2002). Searches for diffuse gamma rays from dark matter annihilation must 


of course also be differentiated from gamma rays that are produced from cosmic rays and 
other astrophysical sources such as AGN and supernova remnants. Both the uncertainties 
in the intrinsic emission and the uncertainties in the astrophysical gamma-ray backgrounds 
provide a substantial challenge for signal detection. 
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III. FERMI-L AT AND I ACT ANALYSES. 

Over the course of the past decade, the field of indirect dark matter searches with gamma 
rays has matured substantially. This is in large part due to the performance of the present- 
generation space and ground instruments, notably the Fermi gamma ray space observatory 
and several lACT telescopes, such as H.E.S.S., VERITAS, and MAGIC. The results from 
these experiments are important both when viewed independently and as results comple¬ 
mentary to those from collider and direct dark matter searches. In this section, after a brief 
presentation of these instruments, we discuss the analysis challenges that they face. 


A. Current gamma-ray instruments 


As the atmosphere is opaque to gamma rays, space borne telescopes are a priori necessary 
to observe the sky at high photon energy. Above a few tens of MeV, their interactions are 
completely dominated by pair production, so that all recent gamma-Oray space instruments 
rely on a pair converter associated to a tracker to detect the trajectory of the produced 
electron and positron pair. The tracker is supplemented with an electromagnetic calorimeter 
to contain the shower and allow for total energy estimate, and an instrumented shield to 
veto the much more frequent incident charged cosmic rays. The most important of such 


instruments currently active is the Large Area Telescope (LAT, Atwood et al., 2009) onboard 


the Fermi space satellite, which also includes a gamma-ray burst detector (GBM). Since June 
2008, the LAT surveys the whole sky every 90 minutes, in the energy range between 20 MeV 
to greater than 300 GeV. The upper end of this range is here determined only by sparse 
statistics for conventional (power-law) astrophysical sources and the availability of reliable 
instrument response functions, i.e. can be much higher. 

The LAT is likely planned for operation until at least 2018. 

Above a few hundered GeV, the fast decreasing flux of typical astrophysical sources and 
the necessarily limited effective area of a space instrument combine to significantly degrade 
the sensitivity. Despite the atmospheric opacity, ground instruments have proven to be able 
to “take over” in this very-high energy range. A similar argument holds for a hypothetical 
WIMP signal, as the flux from WIMP annihilation is inversely proportional to the WIMP 
mass squared. To detect photons at such high energies. Imaging Air Gherenkov Telescopes 
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detect the Cherenkov radiation produced by the charged particle cascade that a high energy 
gamma ray initiates in the upper layers of the atmosphere. The Cherenkov light is reflected 
by large (~ 10 m diameter) mirrors onto cameras consisting of arrays of photomultipliers. As 
the interaction takes place high in the atmosphere (~ 10 km), and the shower needs to reach 
a certain size to be detectable, the analysis threshold is usually close to 100 GeV, though 
larger telescope mirror size (and the corresponding possibility to detect fainter Cherenkov 
light) can push the threshold down to about 30 GeV. The Gherenkov telescopes currently 
in operations are H.E.S.S (|Hinton 2004, 5 telescopes, including the recently added 28- 


meter telescope), VERITAS (Holder et al.. 2006, 4 telescopes) and MAGIG (Lorenz, 2004 


2 telescopes). At slightly higher threshold energy, a different detection technique consists 
in building an array of water tanks instrumented with photomultipliers, and in recording 
the passage of the shower particles themselves, throught the Gherenkov light that they emit 


while traversing the water. This technique is exemplihed by MILAGRO (Atkins et al ., 2003), 
and more recently HAWG (see Section 0. 

The LAT and lAGTs have very distinctive operation modes and face quite different 
challenges. While the LAT surveys the whole sky continually without any need to point in a 
specihc direction (but for the case of transients), lAGTs operate in a pointed mode and with 
a limited duty-cycle (though this is not true for water Gherenkov instruments) and have to 
rely on the larger effective area (by a factor of 10^ to 10^) as compared to the LAT. As a 
result of the held of view, the LAT can a priori provide observations on any of the targets 
discussed in Sectionj^ while in the case of lAGTs, these targets are in competition with other 
astrophysical sources of interest to the community. Furthermore, the LAT has an excellent 
cosmic-ray rejection power, so that its background in conventional analyses is dominated 
by gamma-ray events from various sources in the held of view. The complexity of the LAT 
analysis is due to the multi-dimensional response functions and the need to properly model 
the sky in the region of interest. Less efficient background rejection implies that lAGTs 
are dominated by an isotropically-distributed charged cosmic ray background. In addition, 
the fact that the atmosphere is used as a detector volume implies that varying atmospheric 
conditions as well as the choice of atmospheric conditions and night-sky background induce 
important and difficult to handle systematics in acceptance corrections, which eventually 
(together with the very hard to distinguish component of cosmic ray electrons) provide the 
fundamental limitation to detecting dark matter. 
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B. Astrophysical challenges 


The astrophysical challenges outlined in Section|^call for advanced analysis methods that 
allow robust inference about the particle physics properties of dark matter. Fermz-LAT anal¬ 
yses are largely based on a maximum likelihood method, where spatial and spectral models 
of both signal and background components are fit to data after convolution with a parame¬ 
terized instrument response. For a dedicated review of statistical aspects of these analyses, 
see Conrad (2014). Modeling the background can in certain cases introduce significant 


systematics, with likely the most problematic case being modeling the diffuse gamma-ray 
background towards the Galactic center. Until recently, lACT searches for WIMP dark mat¬ 
ter focused on comparing integral flux predictions with the data, instead of using the full 
spectral information predicted by a single or combination of annihilation or decay channels. 
Modeling of the background has not been commonly performed but rather the background 
expectation was obtained from an off-source region of interest. The signihcance of an excess 


is then commonly inferred from the maximum likelihood ratio test statistic (Li and Ma 


1983) applied to the case of OFF estimation of the background, and confidence intervals 


using either a counting experiment prohle likelihood 

Rolke, Lopez, and Gonrad 

2005 

) or 

using a Neyman construction ( 

Feldman and Gousins 

1998 

). In the remainder of this 

sec- 


tion, we will discuss in detail the different analysis approaches for the different astrophysical 
targets. 

C. Galactic Center 


As discussed in Section the Galactic Genter (GG) is probably the strongest source of 
gamma-ray radiation due to dark matter annihilation. However, the GG is crowded with 
conventional gamma-ray sources: H.E.S.S.and FermULAT sources at the GG are consistent 
with each other and known sources (Gohen-Tanugi, 2009 Gohen-Tanugi et al.\ 2009 Hooper 


and Linden, 2011). For this reason, the GG is typically not directly studied in search for 


dark matter annihilation. Modeling the diffuse emission in the inner Galaxy, defined here 
as the region interior to the Solar circle, also poses a significant theoretical challenge- for 
this reason the Fermi-LAT has not published dark matter constraints from either the GG 
or inner Galaxy, leaving aside searches for spectral features (see below). 
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At lACT energies the systematics of the diffuse emission are less signihcant. The analysis 


presented in Abramowski et al. (2011) applies an ON-OFF technique, where the background 


is determined by OFF-source observations, and in this specihc case the OFF-source region is 
dehned within the held of view. The disadvantage of dehning the OFF-source region within 
the held of view (i.e. relatively close to the GO) is that a potential signal will be subtracted 
as part of the background. This means that the sensitivity might be reduced for relatively 
cored dark matter prohles. An alternative is therefore to dehne the oh region from separate 
pointings, i.e. pointings that are truly oh-source. The likelihood for this approach is: 


noFF\s, b, a) = Pois{noN\s + b)Pois{noFF\c(b) 


( 8 ) 


where uqn and hqff denote the number of counts in ON and OFF regions, respectively, 
a is the ratio between the total ON region acceptance and OFF region acceptance and b 
as usual denotes the background expectation. There are two assumptions going into this 
approach: (a) the background expectation is the same in the on source and oh source 
region, and the ratio between the on-source and oh-source acceptance is known. The latter 
is a systematic which will become important for future lACTs, such as CTA, where the 
statistical uncertainties are sub-dominant as compared to systematic. For a potential way 


of handling these uncertainties see Dickinson and Conrad (2013); Spengler (2015). 


To circumvent the aforementioned disadvantage for the GC that in standard observation 
mode the OFF regions are within the few-degree held of view of the camera, a technique 
which obtains OFF data from a truly separate pointing has to be implemented. However, 
this introduces new systematic uncertainties that have to be carefully addressed; using this 


technique competitive limits have nonetheless been presented in (Abramowski et al. ^ 2015). 


D. Galactic diffuse emission 


The Galactic dihuse emission provides a potentially powerful target for dark matter 
searches, e.g. (Fargion et al., 2000) (Zeldovich et al. , [1980 ), For these searches, the gamma ray 
spectrum and spatial distribution can both be used as discriminants. However, astrophysical- 
induced dihuse emission from the Galaxy is very difficult to model. Parameters to be 
considered are halo height, dihusion coefficients and indices, Alvfen velocity, power-law 
indices of the injection spectrum of cosmic rays and spatial distribution of cosmic ray sources 
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as well as maps of the interstellar radiation field and gas distributions. The most detailed 
model of Galactic diffuse emission and cosmic-ray propagation is provided by the GALPROP 


code (Strong, Moskalenko, and Ptuskin, 2007)^. The more than 20 parameters entering this 
modeling (and in principle also versions of the input gas maps and interstellar radiation 
helds) constitute a set of nuisance parameters that must be handled in order to extract 
a potential dark matter contribution. The most advanced attempt at dealing with these 


parameters is presented in Ackermann et al. (2012a), which uses a prohle likelihood. A 
subset of the full parameter space was considered in this work. For the linear parameters 
a £t to the data provides a prohle estimate, while for the non-linear parameters a grid of 
likelihood points was considered to map the complete likelihood function. However, the 
performance of this mapping is still difficult to assess. The complexity of the problem 
suggests that scanning algorithms designed for Bayesian inference could be more suitable 
for the problem. This has been applied to cosmic-ray data to provide constraints on diffusion 


parameters (Trotta et al, 2011). 


E. Dwarf spheroidals 

For reasons highlighted in Section dSphs are valuable targets for dark matter searches. 
There is now a substantial body of literature devoted to understanding how to extract a 
gamma-ray signal from dSphs. For this reason, we will give here detailed account of the 
methods that are used in these searches. 


1. Formalism and illustration 


Here we describe in detail the full likelihood analysis method used in dSph searches 


with the Fermi-LAT and more recently also with H.E.S.S. (Abramowski et al, 2014b), and 
MAGIG (Aleksic, Rico, and Martinez, 2012) In order to illustrate the advantages of the 


full likelihood analysis, we start by hxing our notation. We consider a dataset Di selected 
in a region of observation “1”, and a statistical model which contains a set of parameters 
of interest p* and a set of nuisance parameters p". We further write the logarithm of 


^ http://galprop.stanford.edu 

^ for an alternative statistical approach see (Geringer-Sameth, Koushiappas, and Walker 20141. 
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the likelihood function £i(Zi)i|p‘, p”). While the maximum likelihood inference based on 
the observation Di results from minimizing Ci with respect to the parameters, it can be 
generalized in a straightforward manner to the case where different, disjoint, observations 
Di and D 2 are modeled using the same parameters of interest. 

The combined log-likelihood function to be used for inference now simply reads 

>Ccomb(^l,^2|p',P?,P2) = >ClPl|p‘,P?) + >C2(^2 |p‘,P2)- (9) 

Such a property of “additivity” is an extremely appealing feature of the likelihood inference 
process, that we now illustrate with the case of a WIMP gamma ray annihilation analysis 
in the direction of two dwarf spheroidal galaxies, each in either region “1” or “2”. The 
parameters of interest are (nn), ttiwimp, and we single out, for the sake of the discussion, the 
J-factors Ji and J 2 from the rest of the nuisance parameters.^ Taking note of the fact that 
the likelihood model in such an analysis actually includes only the product of (crv) and Ji, 2 , 
the combined log-likelihood now reads 

Ccomh{Di, D 2 \{(^v) X Ji,{av) X J 2 , Ji, J 2 ,Pi,P 2 ) = Ci{Di\{av) x Ji, Ji,p") (10) 

-F /:2(^2|(o-t') X J2,-^2,P2)- 

In order to infer a best-£t value for {av) it is thus necessary to fix Ji ^2 to some values; this 
bars any propagation of the corresponding uncertainties into the analysis. This can actually 
be circumvented by considering that any value of a J-factor comes from a set of stellar 
observations d and an inference process based on a posterior distribntion or a likelihood 
function, so that we generically write Pj{d\J,I). Here I stands for any other information 
needed to construct this function, for instance the choice of a parameterized model for the 
dark matter profile. The combined likelihood analysis can be extended to read: 

Ccomh{Di,D 2 ,di,d 2 \{(Jv) X Ji,{av) x J 2 , Ji, J 2 , p", P 2 ) = Ci{Di\{av) x Ji,Ji,p^) (11) 

+C2{D2\{av) X J2,^2,P2) 

+Pj{di\Ji, I) + Pj{d2\J2,1)- 

As a result the degeneracy between the J-factor and (av) is removed, and there is an auto¬ 
matic propagation of the statistical uncertainty in the stellar analysis down to the inference 

^ In practice, the WIMP mass is kept fixed during the likelihood inference process, so that we omit it in 
the rest of this discussion. 
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about {(Tv). Given the shape of the Bayesian posterior probability densities commonly ob¬ 


tained in J-factor derivations from stellar analyses (Martinez et al., 2009 Strigari et al. 


2008), Ackermann et al. (2011[) actually use a Log-Parabola function as an ansatz for the 


Pj function. 

We can now illustrate the power of this combined analysis. For this purpose we consider 
the two dSphs, Carina and Ursa Major II. For those two dSphs, the central and sigma 
values for the Log-Parabola function are (6.3 x 10^^, lO*^'^) and (5.8 x 10^^, 10°'^), respectively. 
Thus Carina has a much lower J-factor than Ursa Major II, but its distribution is much 
narrower, so the uncertainty incurred in choosing the central value as hducial J-factor is 
much lower. Figure presents the likelihood curves for the single-dSph and the combined 
analyses, prohled over the nuisance parameters, including the J-factors. 



(a v) ( xlO ) 


FIG. 2 Relative likelihood function versus velocity averaged annihilation cross-section. The effect 
of a combined likelihood analysis with inclusion of the J-factor uncertainty, in the case of two 
dSphs, Carina and Ursa Major II. is illustrated. 


The mass and channel used for this illustration are 100 GeV and bb, respectively. The 
intersection between the horizontal dotted curve and a likelihood curve localizes in the x- 























axis the 95% C.M. value for (cru), so that the lower its value, the stronger the upper limit. 
Thus, as expected when no account is taken of the J-factor uncertainties, Ursa Major II 
(green dashed curve) outperforms Carina (blue dashed curve) by and large. On the other 
hand, the combined analysis (red dashed curve) performs barely better than Ursa Major II 
alone, due to the modest statistical power of Carina. When J-factor uncertainties are taken 
into account (see plain curves, with identical color coding), one immediately notices that 
Carina’s limit does not degrade by much, while limits from Ursa Major II are signihcantly 
weaker than before. And furthermore, the combined likelihood result is now signihcantly 
improved compared to Ursa Major II ; when J-factor uncertainties are taken into account, 
the statistical power of Carina is no longer negligible. Comparison of the dashed and plain 
red curves shows the effect of the J-factor uncertainty on the combined upper limit. 


2. Dwarf spheroidals and the global fit 


The likelihood method described above can be extended to attempt statistical inference on 


a specihc supersymmetric model, including taking into account collider experiments (Scott 


et al, 2010). The result of this type of analysis is a likelihood or posterior maps in the 


parameter space of constrained supersymmetry. However, because of the sparsity of the 
data, inference on the physical model is difficult. Nonetheless, this result serves as a proof of 
concept for future attempts to combine results of different dark matter probes in a consistent 
likelihood or Bayesian inference. 


F. Dark satellites 

Dark satellites are intriguing objects for dark matter searches. In an ideal scenario, they 
would constitute a gamma ray source without counterpart in any other wavelength. For a 
preliminary estimate, the list of unassociated Fermz-LAT sources provides a catalogue of 
potential dark matter satellites. 

In order to be classihed as a dark satellite candidate, hrst and foremost the energy 
spectrum should be consistent with a dark matter induced spectrum and distinguishable 
from more common astrophysical-induced spectra, such as a power-law, or a more difficult 


to deal with pulsar spectrum (see e.g. Baltz, Taylor, and Wai, 2007). Additionally, the 
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source may be extended, which would be in particular true if the source were very near to 
the solar system. From a statistical perspective, a likelihood ratio test statistic provides a 
mean to address this question of whether the energy spectrum of the source is consistent 
with dark matter. The problem though is that the hypotheses, in particular the question 
whether the source exhibits a power-law as compared to a dark matter spectrum, constitutes 


a non-nested model comparison for which the usual asymptotic theorems (such as Wilks 


1938) do not hold. This implies that the signihcance of a potential detection (rejection 


of the null hypothesis) cannot be robustly calculated. Techniques that are proposed to 


address this problem are known in literature (James, 2006), but it is unclear if they perform 


in a satisfactory manner (see e.g. Conrad, 2014, and references therein). Thus, the null 


distribution of the test statistic is derived from Monte Carlo simulations of the experiment 


or from random region of interests in the sky (Ackermann et al, 2012c). For lACTs, the 


list of unassociated sources provides potential targets for deep follow-up observations, and 
the same criteria for claiming a detection would apply. 

G. Galaxy Clusters 


The search for dark matter in Galaxy clusters utilizes analysis methods similar to those 
used for dSphs. The main difference is that galaxy clusters are in all likelihood extended 
sources of gamma ray emission. Indeed, if substructures and their corresponding flux dom¬ 
inates the dark matter contribution, it is the outer parts of clusters that dominate. As a 
specihc example, the Virgo cluster may have an extension of approximately 6 degrees. The 
search for gamma ray emission from clusters therefore has to use extension as source model. 
The major drawback of this is that the results are in general sensitive to the modeling of 
the diffuse emission, and uncertainties in this diffuse emission must be accounted for. A full 
mapping of the likelihood space of this component is very difficult (see also next subsection). 


As a result, in the analysis presented in Ackermann et al. (2014a) the systematic effects due 


to modeling of the diffuse emission are accounted for by considering a set of fiducial diffuse 
models and recalculating limits separately for each of these fiducial models. 
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H. Searches for dark matter line signal 


Searching for spectral features is in principle a simpler task than searching for a continuum 
signal. In this case, the main question is whether there is a local (in energy space) deviation 
from a background, which can be inferred from the data itself. This means that no physics 
modeling of the background is necessary. The line signal is described by a delta function 
convoluted with the energy dispersion introduced by the instrument. 

The analysis of line signals typically proceeds by applying a sliding window technique, i.e. 
the window is centered on the particle mass to be tested, and the background is determined 
either from the whole window or excluding some signal region. The size of the window 
is chosen such that the background is described by a simple empirical £t (in best case a 
power-law). Statistical inference is often performed using an extended likelihood £t, with 
application of Wilks theorem ( Wilksf 1938) or Chernoff theorem (Chernoff, [1954 ) to obtain 
signihcance and the prohle likelihood to infer conhdence intervals. At this stage, most 
analyses do not attempt to model the spatial distribution of the gamma rays, but instead 
optimize the region of interest using a signal and background prediction and then treating 
the problem only with the spectral likelihood. The above approaches have been applied 
in 


Abdo et al. 

(2010a) 

Ackermann et al. 

(2012b 

2013) 

Weniger 

(2012) 

Ackermann et al. 


(2013) also takes into account the quality of the event reconstruction by considering not 
only the likelihood of the parameter of interest but by also considering the distribution of 


photons in a quality parameter. The analysis in Su and Finkbeiner (2012) uses both spectral 
and spatial distributions. 

While peak hnding of this sort is as old as particle physics, challenges remain. Potential 
challenges come from choosing the side band window large enough to allow a good determi¬ 
nation of the background, while still being small enough that a simple empirical function will 
describe it. Another important aspect is that very good knowledge of the energy dispersion 
is necessary as inference is done on a steep spectrum where small biases in our knowledge 
of the energy dispersion can have large effects. Finally, in searches for dark matter, mass 
is a parameter that is not dehned under the null hypothesis and will thus introduce a trial 
correction, which in practical applications is not simply calculated by applying the binomial 
distribution, but correlations have to be taken into account. In these more complicated 
cases there are two approaches to follow: employing Monte Carlo simulations or resampling 
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from off-source data distributions. The main challenge here is to be able to simulate the null 
distribution to sufficient accuracy for the high signihcance needed (usually 5 a), i.e. the need 
for a large (~ 10®) number of independent experiments. Solutions to this have emerged, the 
simplest being to £t an empirical function to the obtained test statistic distributions (e.g. 


done by Weniger, 2012). Another one proposed by (Gross and Vitells, 2010) extrapolates to 


higher significance from a lower number of simulations (see also Conrad, 2014, for a more 
detailed discussion). It should be noted that the latter is only strictly applicable in the case 
of continuous trials, i.e. the whole spectrum is fit with the mass as a free parameter. Finally, 
in case of not very conclusive data, it might be difficult again to distinguish a line feature 
from any other feature (e.g. a broken power-law). Again, non-nested hypothesis testing 
would have to be performed for rigorous results. 


IV. STATUS OF DARK MATTER SEARCHES WITH GAMMA RAYS 


In this section we review the current status of indirect dark matter searches with gamma 
rays. We assimilate the information on astrophysical targets (Section 0 with the infor¬ 
mation on instrumental sensitivities and statistical methodology (Section |III| ). For the two 
main types of gamma ray experiments, the all-sky Fermi-LAT and the lACTs, our goal is to 
review the statistical analysis method that is most appropriate for each target, and review 
the results that have been obtained. 


As was reviewed in Section the nominal all-sky survey mode of the Fermi-LAT is 
ideally suited to explore essentially all potential targets in the gamma-ray sky. From an 
astrophysical and experimental perspective, though, these targets incur signihcantly different 
systematic uncertainties. 

As we detail in this section, in the past few years the Fermi-LAT collaboration has 
extended the statistical framework of the official ScienceTools to include a joint likelihood 
formalism that allows for the combination of several source regions of interest into a single 


analysis. Such a combined analysis was introduced in Ackermann et al. (2011), and was 


reviewed in Section III.E above. Though this framework is not yet able to completely solve 
the complexity of combining LAT analyses of all possible targets, it has for the hrst time 
produced sensitivities to the nominal thermal relic value of {av) ~ 3 X 10 cm®s 

At higher energies, lACTs have gained ground. The constraints obtained from observa- 
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tions of the Milky Way halo in the vicinity of the Galactic Center constitute a large step 
forward for constraining WIMP masses above about a TeV. 

In the following sections, we provide a summary of the most relevant gamma ray analyses 
and results from Fermi LAT and lACTs, including observations of the Galactic center 
and diffuse Galactic halo, dSphs, dark satellites, galaxy clusters, and the diffuse isotropic 
background. 

The spectrum of gamma rays depends on the annihilation channel. Unless otherwise noted 
we will follow the convention to present results for annihilation into b-quark pairs, which 
can be seen as representative for quark annihilation in general. Annihilation to r-leptons 
generically results in a harder spectrum. Majorana WIMPs in general do not preferably 
annihilate to light leptons (electrons or muons) as the annihilation is helicity suppressed, 
except if for example virtual internal bremsstrahlung or hnal state radiation is considered 


for supersymmetric particles, e.g. (Bringmann, Bergstrom, and Edsjo, 2008). 


A. Galactic Center 


As a prime target for indirect dark matter searches, the Galactic center has been the 


subject of intense scrutiny for well over a decade. Gesarini et al. (2004) proposed a dark 


matter interpretation to the gamma ray source detected with EGRET (Mayer-Hasselwander 


et al, 1998). Near the same time lAGTs also detected a TeV source (Aharonian et al 


2004 Albert et al, 2006), which seriously challenged a dark matter interpretation of the 


EGRET source (Aharonian et al, 2006). Independent of dark matter, these observations 


were important because they showed that near the Galactic center, in addition to the diffuse 
background, there exist strong gamma ray emitters. More recent analyses of the Fermi- 


LAT data (e.g. 

Abazajian et al, 

2014alb; 

Abazajian and Kaplinghat, 

2013; 

Galore, Gholis, 

and Weniger, 2014; Daylan et al. 

2014 

Hooper and Goodenough 

2011 

) have reported the 


detection of an extended excess, compatible in shape and spectrum with a ^ 30 GeV WIMP 
annihilating to bb. 

Several groups have followed up and conhrmed the existence of this inner Galaxy excess 


in the Fermi-LAT data. Abazajian et al. (2014a) pointed out the critical dependence on 


foreground modeling, as well as the possibility that the excess arises from unresolved mil¬ 


lisecond pulsars. Garlson and Profumo (2014); Petrovic, Serpico, and Zaharijas (2014) show 
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that this excess can be explained by local population of cosmic-ray protons, potentially due 
to a burst-like injection several thousand years ago. There is also the yet unsettled char¬ 


acterization of the Fermi bubbles at low latitude (see Galore, Cholis, and Weniger (2014) 
for a discussion of this point). At the time of writing, there is no clear resolution to this 
excess emission, and it is quite possible that there won’t be a resolution for some time. It 
is certainly true that a conhrmation of this excess in one of the other targets, for example 
clusters or dwarf galaxies, would be needed to further strengthen the case for dark matter. 
Taking a best guess nominal J-factor of the GC, the most recent Fermi-LAT combined 


dSph search could have conhrmed the excess, but only reported constraints (Ackermann 


et al, 2015) (see below). Probably the most conservative stance at this point assumes that 


all of the excess photons are not due to dark matter annihilation; this assumption implies 


stringent constraints on a contracted NFW scenario at the Galactic Genter (Gomez-Vargas 


et al ., 2013). 


B. Galactic Halo 


The diffuse Galactic halo has long been advocated as an interesting target for dark mat¬ 


ter searches. Serpico and Zaharijas (2008); Stoehr et al. (2003) hrst discussed the idea of 


searching in an annulus around the center of the Galaxy. This would bypass aforementioned 
systematic uncertainties that hinder such a Galactic center analysis. The HESS collabora¬ 


tion put such a strategy in practice in Abramowski et al. (2011), where an annular region 


from 0.3° to 1° radius was analyzed, excluding the Galactic plane (|6| < 0.3°). Using 112 


hours of livetime observations and a generic continuum spectrum from Tasitsiomi and Olinto 


(2002), they derived, for NFW or Einasto density prohles, the most stringent limits to date 


in the 1 TeV WIMP mass range. For instance (av) in excess of 3 x 10 cm^s ^ is excluded 
assuming an Einasto density prohle. 


As discussed in Section dehning the background region within the held of view is not 
efficient, especially if the dark matter distribution in the central galaxy is cored. A proof 


of concept for treating this situation is published in (Abramowski et al, 2015) providing 


competitive limits with only 9 hours of data. The main conceptual novelty is that the 
OFF region is constructed from separate pointings of the telescope array. To make sure 
that atmospheric conditions are preserved as much as possible an offset in right ascension 
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is chosen to define the OFF-source pointings. However even in this case a novel treatment 
of the acceptance correction for the background estimate had to be introduced. In this 
particular case the center of the OFF regions is about 6 degrees away from the center of the 
signal region, and the J-factor is reduced by at least a factor of 3. 

More recently, Fermz-LAT has analyzed data on a larger scale within a region of |/| < 15° 


and 5 < l&l < 80° (Ackermann et al., 2012a). This analysis used Pass 7 and masked 
IFGL point sources. The analysis proceeds with the CLEAN class and thus the P7CLEAN_V6 
instrument response functions (irfs) to infer stringent limits on (crv), while discussing in 
great details the dependence of the analysis on the Galactic foreground modeling, which is 
the main issue with such a halo analysis when comparing to the dwarf spheroidal constraints. 

To summarize the previous discussion and the current observational situation, we gather 
in Figure [^the most relevant current limits obtained for analyses of regions within the Milky 
Way. This figure clearly illustrates the power that a deep observation with a Gherenkov 
instrument has to set competitive upper limits for energies about ~1 TeV. It also emphasizes 
the advantage of a halo analysis with the LAT, with respect to satellite or GG blind searches, 
though we stress again that background systematics affect the limits much worse in this case. 
In addition, the halo upper limits are proportional to the squared normalization of the dark 
matter density distribution of the Galaxy, usually estimated at the position of the Sun. As 
discussed in Section there is still a substantial uncertainty in the estimation of the local 
dark matter density, 0.2 to 0.9 GeV cm“^, so that the resulting uncertainty on a dark matter 
annihilation flux is in fact larger than the uncertainty just contributed by the modeling of 


the Galactic diffuse emission (see Ackermann et al., 2012c, for further details). 


C. Dwarf spheroidals 

As discussed in Section |TT[ with their lack of high-energy emission processes and large 
dark to luminous mass ratio, dwarf spheroidal galaxies are a prime target for gamma ray 
instruments. As the expected flux scales with the inverse square of the distance, only the 
satellites of the Milky Way have been investigated so far, which amounts to about 25 targets. 
While Fermi is surveying the whole sky continuously, lAGTs need to allocate observation 
time to point in the direction of any dSph. As of this work, a total of 9 dSphs have been 
observed by AGTs. The HESS collaboration initially presented constraints from ~11 hours 
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FIG. 3 Upper limits (UL) on velocity averaged annihilation cross-section versus WIMP mass from 
indirect searches in 7 rays, for analyses focusing on the Milky Way. In dashed red (“HESS halo”), 


the 95% CL UL in the Galactic Center halo analysis from Abramowski et al. (2011); in red (“LAT 


halo no bkg”), the LAT halo analysis (Ackermann et a/., 2012) corresponding to 3cr CL upper 
limits, when no diffuse background modeling is performed; in green (“LAT halo”), same but in the 
case where proper modeling of the diffuse gamma ray background is performed; in cyan (“LAT 


GC”), the LAT diffuse-model-free 3 cj ULs at the Galactic Center (Gomez-Vargas et al.. 2013); 
in magenta (“LAT GC rescaled”), same but rescaled to the same local DM density as the LAT 
halo analysis; in yellow (“LAT GC NFWc), same but in the case of a contracted NFW profile; in 


magenta (“LAT HVC”), the 95% C.L. limits obtained with the Smith cloud (Drlica-Wagner et al. 


2014). The dot marker with the label “satellites” corresponds to the 95% CL upper limit for a 100 


GeV WIMP mass obtained in the unidentified LAT source analysis (Ackermann et al.. 2012c). 


of observation toward Sagittarius dwarf (Aharonian et al, 2008, 2010), which have recently 
been largely revised with 90 hours accumulated on target and with a major modihcation 


to the density prohle of this tidally stripped nearby galaxy (Lamanna et al.. 2013). The 
HESS collaboration also published results from Canis Major (~ 10 hours) (Aharonian et al 


2009c), Sculptor (~ 12 hours) and Carina (~ 15 hours) (|Abramowski et 'd., 2011). As 


mentioned previously, HESS also has presented an analysis combining previous observations 
of dSphs with a new long exposure of the Sagittarius dwarf, employing a combined likelihood 


technique inspired by the Fermi-LAT analysis (Abramowski et al.. 2014b). The MAGIC 


collaboration presented results from 10 to 20 hours of observations towards Willman I (Aliu 



















































et al ., 2009), and Draco (Albert et al ., 2008), but devoted most of the allocated time to the 


observation of Segue 1, initially with 43 hours on one telescope (Aleksic et al, 2011), and 


recently updated with a 160-hour analysis (Aleksic et al, 2014), that also makes use of a full 
likelihood technique, akin to what the Fermi-LAT collaboration is routinely doing. Finally, 
the Veritas collaboration presented results for Bootes I, Draco, Ursa Minor and Willman 


I (Acciari et al, 2010), and for Segue 1 (50 hours) (Alin et al, 2012). 

These observations are to be compared to the 1000-hour equivalent of 11 months of 
LAT survey, which resulted in flux limits set in the direction of 20 dSphs, and dark matter 


constraints derived for a subset of 8 of these, based on robust J-factor values (Abdo et al 


2010c). Focusing specihcally on the constrained minimal super-symmetric model, Scott et al 


(2010) added the analysis of Segue 1 with approximately the same amount of data. Thus, 
quite generically given the plausible allocation time on dSphs by lACTs, LAT constraints 
are expected to dominate below WIMP masses of order 1 TeV, at which point the fast 
degradation of LAT sensitivity quickly limits its performance. As a consequence, while 
programs are ongoing to increase the total duration of lACT observations in the direction 


of dSphs (e.g. Zitzer et al, 2013), it is also crucial to maximize the statistical power of the 
analyses by using the full-likelihood technique, as considered for instance in Aleksic et al 


(2014); Aleksic, Rico, and Martinez (2012) and Abramowski et al (2014b). 


Along these lines. Acker ma un et al (2011) introduced a joint-likelihood formalism that 
combines the statistical power provided by several dSphs into a single inference (see Sec¬ 


tion III.E, and also Geringer-Sameth and Koushiappas (2011) for an alternative combined 
methodology). Using 24 months of Pass 6 “diffuse” data and the corresponding P6V3 irfs, 
the Fermi-LAT collaboration derived a single 95% C.L. exclusion curve using 10 dwarfs, 
that for the hrst time reached the “standard” thermal-relic value of (av) = 3 x 10“^® cm“^ 
s“^, for annihilation into quarks. Furthermore, the Fermi-LAT collaboration introduced for 
the hrst time a scheme to directly account for the J-factor uncertainty into the upper limits 
on (av) stemming from the statistical nature of the dark matter prohle derived from stellar 


data. The analysis presented in Ackermann et al (2011) was later updated to 4 years of 


Pass 7 reprocessed data (Ackermann et al, 2014b), in which a much more thorough analysis 
of the data in the direction of 25 dSphs was done. A subsample of 15 dSphs with robust 
J-factors were retained to derive upper limits on the annihilation cross-sections in different 
channels. In the latter analysis, it was found that the test statistic distribution, naturally 
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defined as the likelihood ratio between the null hypothesis (background-only) best £t and 
the DM-|-background £t did not follow a distribution, implying a correction reducing the 
apparent signihcance of an excess. The most likely cause of this deviation is a population 
of unresolved point sources, confirmed by the most recent incarnation of the Fermz-LAT 


analysis (Ackermann et ai, 2015). 


For the case of dSph searches it is worth noting that the advantage of the full likelihood 
formalism developed lies in the fact that it opens up the possibility to eventually share the 
likelihood functions across collaborations to add as many dSphs as possible, or even com¬ 
bine targets. As a hrst step in this direction, the Fermi-LAT collaboration has released the 


likelihood functions used in Ackermann et al. (2014b). Finally, to our knowledge the possi¬ 


bility to add J-factor uncertainties to a full-likelihood formalism, as proposed by the LAT 


collaboration (Ackermann et al.. 2014b Aleksic, Rico, and Martinez, 2012), has only been 


considered in the most recent search for dSph emission performed by the HESS collaboration 


(Abramowski et al, 2014b) 


On the topic of satellite galaxies of the Milky Way, it also should be noted that the 
Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC) also provide targets for 
dark matter searches. Extraction of gamma-rays from dark matter must however compete 


with the detected gamma-ray emission that arises from cosmic rays (Abdo et al, 2010d 


Abdo et al., 2010c). Buckley et al. (2015) have recently performed an analysis of the LMC 


searching for a dark matter signal with FermALAT data, and report a null detection. 


D. Dark satellites 


The FermALAT has performed a search for a gamma-ray signal from the predicted pop¬ 
ulation of dark satellites of the Milky Way ( Ackermann et a7(| 2012). In this analysis the 
231 high-latitude (|5| > 20°) unidentihed IFGL sources were augmented with 154 candidate 
detections using a dedicated search for potentially extended sources at high latitude, and 
one year of Pass 6 DIFFUSE-class LAT data (in the energy range 200 MeV to 300 GeV). A 
spectral and spatial selection was applied to check compatibility with a dark matter signal, 
resulting in no remaining candidate and, after comparison with simulations, a constraint on 
{(Tv) is derived which is about 1.95 x 10“^^ cm^s“^ for a 100 GeV WIMP annihilating into 


hh. Starting from the IFGL point source catalog, Buckley and Hooper (2010) also undertake 
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a search for dark matter subhalos, and similarly report no conclusive detection. 

Some theories suggest that high-velocity clouds (HVC) may be embedded within dark 


matter subhalos (Blitz et a/., 1999), and if so, they provide a target for gamma ray 


searches. Drlica-Wagner et al. (2014) focused on the Smith Cloud, a massive low-metallicity 


HVC located at a distance of about 12 kpc from the Sun (this is one of the few HVCs with a 
known distance). Five years of Pass 7 reprocessed CLEAN-class data in the energy range 500 
MeV to 300 GeV yielded no detection, resulting in constraints on (av) which are compara¬ 


ble to the dSph combined analysis (Ackermann et al., 2014b) for an assumed NFW prohle. 


Of course it should be emphasized that the dark matter content of HVCs is very contro¬ 


versial (Faerman, Sternberg, and McKee, 2013), and the constraints from the Smith Cloud 


strongly depend on the dark matter prohle. Indeed the limits degrade by a factor 40 when 


using a Burkert prohle (Burkert, 1995) rather than a NFW or Einasto prohle (Drlica-Wagner 


et al ., 2014). 


E. M31 


In much the same way as the Milky Way, the Andromeda or M31 galaxy, our closest spiral 
galaxy neighbor, is expected to shine in gamma rays. The primarily emission is due to its 
gas and cosmic-ray content; these components provide a background from which a potential 
dark matter signal must be disentangled. In 2010 the Fermz-LAT collaboration announced 
the detection of M31 using about two years of Pass 6 DIFFUSE-class (P6DIFFUSE_V3 irfs) 


data in a 10° x 10° squared region centered on M31 (Abdo et al, 2010a). Using a spatial 


template derived from the IRIS 100 pm far infrared map (Miville-Deschenes and Lagache 


2005), the spectrum derived from a ht to the LAT data is consistent with a rescaled spectrum 


obtained by a GALPROP run to model the Milky Way gamma ray emissivity (Vladimirov 


et al.. 2011). Thus the LAT detection of M31 is compatible with star formation rates and 


gas content, and a conservative 95% G.L. upper limit on {av) for a 100 GeV WIMP in the 
bb channel is derived at abont 5 x 10“^® cm^s“^, assuming a smooth Einasto density prohle 


derived from Klypin, Zhao, and Somerville (2002) 


The M31 analysis has since been updated with a 4.5 year Pass 7 dataset (Li, Yuan, and 


Xu, 2013). These authors add a dark matter prohle to the model that already includes a 


template for the emission from M31 due to its gas content. This makes the resulting upper 
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limit less conservative but the analysis more sensitive. Under the assumption of a smooth 
halo, the cross section upper limit is very similar, ~ 5 x 10“^® cm^s“^. More recently there 


has also been the claim of an extendend excess attributed to a cosmic-ray halo (Pshirkov, 


Vasiliev, and Postnov, 2015). 


Finally, in the TeV energy range, early searches have been performed with CE¬ 


LESTE (Lavalle et al., 2006), and HEGRA (Hofmann et al., 2003), without detection and 


with uncompetitive upper limits derived. 

Figure [^summarizes the status of the dark matter searches with gamma rays in the Local 
Group. Similar to the constraints obtained from the Milky Way, this hgure clearly indicates 
the complementarity between Fermi-LAT and lAGT constraints. 

F. Galaxy clusters and isotropic emission 

1. Galaxy clusters 


As discussed in Section II.E galaxy clusters are anticipated to be gamma ray sources. 


Using EGRET data, no detection was reported in Reimer et al. (2003) for 58 clusters selected 
from an X-ray-bright sample, resulting in an average 95% G.L. flux upper limit of ~ 6 x 10“® 
cm“^ s“^ above 100 MeV. Null detections were also reported above 400 GeV in the direction 
of the Perseus and Abell 2029 clusters with the Whipple telescope, using ~14 and ~6 hours 


on source, respectively (Perkins et al., 2006). More recent null results in searches with AGTs 


include Aharonian et al. ( 2009a||b ); Aleksic et al. (2012, 2010); Arlen et al. (2012). With 


an 18-month generic gamma ray analysis toward 33 clusters (Ackermann et al., 2010a) that 


improved upon previous analyses (Reimer et al., 2003), the LAT collaboration presented 


early constraints on a dark matter induced gamma ray signal in a subset of six clusters 


and for bb and channels (Ackermann et al., 2010b). In the latter case, the cluster 


analysis conhrmed the tension between the LAT e’''e spectrum and a generic leptophilic 


dark matter scenario that would aim to explain the Pamela positron excess (Adriani et al 


2009). The bb limits also showed promise, especially when accounting for expected but 


unknown cluster substructures. Fornax, the best target among the six clusters studied 
in 


Ackermann et al. (2010b), was also observed by HESS with 14.5 hours and a threshold at 


100 GeV (Abramowski et al, 2012, 2014a). The Veritas collaboration also presented dark 
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FIG. 4 Comparison of constraints on the velocity averaged annihilation cross section (bb channel) 
versus WIMP mass derived from the LAT combined analysis of 15 dwarf galaxies (assuming an 
NFW profile), 160-hour observations of Segue 1 by MAGIC (|Aleksic et al. 2014), 48-hour obser¬ 


vations of Segue 1 by VERITAS (assuming an Einasto profile) (Aliu et aZ., 2012), and 112-hour 


observations of the Galactic center by HESS, assuming an Einasto profile (Abramowski et al 


2011). In the interest of a direct comparison, we also show the LAT constraints derived for 


Segue 1 alone assuming an Einasto dark matter profile consistent with that used by VERITAS 


(Aliu et al. 2012). Eor this rescaling, the J-factor of Segue 1 is calculated over the LAT solid 
angle of AH ~ 2.4 x 10““^ sr and yields a rescaled value of 1.7 x 10^® GeV^ cm“® sr (uncertainties on 
the J-factor are neglected for comparison with VERITAS). The Pass 8 limits (in magenta) and 


related expected-sensitivity bands are from Ackermann et al. (2015). Finally, the green dashed 


curve is from the 5-dwarf combined analysis of HESS data by Abramowski et al. (2014b). 


matter constraints in Arlen et al. (2012), based on 18.6 hrs of observations of the Coma 


cluster, also studied in Ackermann et al. (2010b). 


As updated gamma ray constraints have been recently derived by the LAT collaboration 


m 


Ackermann et al. 

(2014a 

), using a combined analysis akin to 

Ackermann et al. 

(2014b), 


it can be expected that a combined dark matter analysis of all or a subset of the clusters 
considered in this paper will soon be presented. For now we can only rely on the preliminary 
combined results obtained with 2 years of Pass 6_V11 DIFFUSE data that were presented 


in Zimmer, Conrad, and Pinzke (2011). In addition to Coma and Fornax, this analysis used 
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M49, Centaurus, and AWM 7. These recent results are illustrated on Figure Finally, 


recent claims of an extended emission toward the Virgo cluster (Han et al.. 2012) have 


not been confirmed (Macias-Ramirez et al, 2012), emphasizing the importance of correctly 


deriving the background model before drawing conclusions based on the Fermi-LAT data. 
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FIG. 5 Constraints on the velocity averaged annihilation cross-section versus WIMP mass from 
analyses in the directions of galaxy clusters, for an NFW profile and bb channel. The red and blue 
curves show the 95% C.L. upper limits obtained with Fornax and Coma clusters, respectively, and 
with Fermi-LAT (thick, Ackermann et al (2010b[)) and lACTs (dashed, Abramowski et al, 2014a 


Arlen et al., 2012|) instruments. The cyan dashed curve shows the 95% C.L. upper limits obtained 


with a combined analysis of 5 clusters (Zimmer, Conrad, and Pinzke, 2011). The cyan thick line 


shows the upper limits derived in Ajello et al. (2015) using the recent LAT EGB observations 


(Ackermann et al.. 2015a). The grey band shows the corresponding uncertainties rising from 


modeling the expected dark matter signal. 


2. Isotropic signal 


Beneath the conspicuous Galactic diffuse emission and the constellation of resolved 
gamma ray sources, which are mostly active galactic nuclei (AGN), the gamma ray sky 


harbors an isotropic signal, the isotropic gamma ray background (IGRB) (Fornasa and 


Sanchez-Gonde, 

2015 

). This was first detected already by EGRET (Sreekumar et al. 

1998), 


and recently determined with very good precision by the FermGLAT collaboration from 


42 


















































approximately 100 MeV - 800 GeV (Ackermann et al ., 2015b). The IGRB certainly includes 


unresolved contributions from standard astrophysical sources, notably AGN. This emission 
may also contain the so-called “cosmological” dark matter signal coming from the summa¬ 
tion of the dark-matter annihilation contributions of all dark matter halos across the history 
of the universe. Approaches to hnd this signal are based on the spectrum of isotropic 


component (see Ullio et al., 2002) or on a spatial signature exploiting the fact that dark 
matter-induced anisotropies in the emission should follow the square of the mass density, 
whereas conventional astrophysical sources should follow the dark matter density linearly. 


which would reveal itself in differences in the angular power spectrum (see e.g. Guoco et al 
2008). Following the hrst measurement of the EGB^ spectrum with Fermi LAT data (Abdo 


et al., 2010b), the LAT collaboration has presented the limits obtained from the spectral 


shape alone (Abdo et al ., 2010b). These constraints are very model dependent in two ways. 


First, the halo and subhalo abundance has to be modeled as function of redshift, and second 
the contribution of conventional sources has to be modeled. In a situation with a poorly 
modeled background a conservative way to place model constraints is to assume that the 
entire detected emission is from the putative signal. With this approach, constraints are 
about two orders of magnitude above the thermal limits even under moderately optimistic 


assumptions on substructure properties (Abdo et al ., 2010b). A proper background modeling 


would improve the constraints by maybe one to two orders of magnitude, see e.g. Bringmann 


et al. (2014). Thus, the publication of an updated IGRB spectrum by the LAT collaboration 


(Ackermann et al., 2015a) comes in par with an upcoming reassessment of the contribution 


of astrophysical sources, notably blazars (Ajello et al, 2015), to the EGB budget. Together 


with an improved modeling of the dark matter expected signal ( |Ackermann et al. . 2015c), 
this effort leads to potentially competitive limits, as shown on Figure]^ 

Finally, anisotropies in the gamma-ray sky have been investigated in [Hensley, Siegal- 
Gaskins, and Pavlidou ([2010 ) and Guoco et al. (2011). Gombining the spatial signal and 


spectral signals may enhance the dark matter sensitivity of these searches (Gampbell and 


Beacom, 2013) and combining isotropic emission measurements with other sources can help 


to break the degeneracy between substructure boost and annihilation cross-section (Ng et al 


2014). 


^ Extra-Galactic Background : the sum of the IGRB and resolved LAT extra-galactic sources. 
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Though employing anisotropy studies in lACTs have potential for very competitive con¬ 


straints (Ripken et al, 2014), so far these analyses techniques have not been applied to our 


knowledge. Systematics from the combination of several fields of view might hamper the 
usefulness of this approach. 


G. Searches for spectral features 

Dark matter lines were suggested as a smoking gun signal almost thirty years ago 


(Bergstrom and Snellman, |1988 ). Experimentally (as pointed on in Section III), line de¬ 
tection can be established relatively independently from background modeling, and source 
confusion is unlikely. Therefore, in contrast to continuum searches, the most promising tar¬ 
get is the Galactic center or its vicinity. Constraints on line emission have been presented 


using EGRET data (Pullen, Chary, and Kamionkowski, 2007)(Mack et a/., 2008), Fermi 


LAT ( 

Abdo et al.\ 

2010a 

Ackermann et al. 

2012b 

Albert et al ., 

2014 

Ger inger- S amet h 

and Koushiappas 

2012 

) and H.E.S.S. data ( 

Abramowski et al.. 2013). A summary of the 


most relevant present upper limits on line emission is shown in Figure 


Weniger (2012) and Su and Finkbeiner (2012) claimed detection of a line emission in 


Fermi-LAT data. Using an analysis technique similar to Abdo et al. (2010a), as described 


in Section III, but with a doubling of the amount of data as well as an optimization of the 


region of interest for signal over square-root of background, Weniger (2012) found a (trial 


corrected) 3.2 a significant excess corresponding to a dark matter mass of ~ 130 GeV. If 
interpreted as a signal, this would amount to a cross-section of about (av) ~ 10“^^ cm^s“^. 
However, this analysis was based on a mis-calibrated data set, and the inferred position 
of the line feature changes to ~ 133 GeV with a correctly-calibrated sample. The claimed 
signal is concentrated on the Galactic center with a spatial distribution consistent with an 


Einasto profile ( 

Bringmann and Weniger 

2012 

). - 

upper limit presented in 

Ackermann et al. 

(2012b) 


In Ackermann et al. (2013), about 3.7 years of Fermi-LAT data were analyzed, including a 


careful study of systematic uncertainties and a discussion of the alleged line feature. Apart 
from using correct calibration, this analysis also employed event-by-event reconstruction 
quality to improve the energy dispersion description. In terms of the line feature this analysis 
also hnds an excess at 133 GeV, but the global signihcance is estimated to be only l.ha when 
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FIG. 6 Upper limits on the velocity averaged annihilation cross section into gamma-ray pairs versus 
WIMP mass, derived by several instruments in searches for line features. Blue and yellow lines 


show the limits obtained with the FermULAT in a first standard analysis (Ackermann et al., 2013) 


and in a dedicated low-energy follow-up (Albert et al., 2014), respectively. The red line shows the 


limits obtained by the HESS collaboration (Abramowski et al. 2013). For reference, the EGRET 


limits are also shown (Pullen, Chary, and Kamionkowski, 2007) 


including information on the event-by-event energy reconstruction. 

Since the hrst report of a line-like feature, the main challenge to the claim that a line 
feature originates from dark matter annihilation was the fact that an excess at the same 
position also appears in gamma rays observed from the Earth’s limb. The signal there 
appears after an appropriate zenith cut is applied that accounts for the fact that the Fermi- 
LAT is usually not ‘head-on’ exposed to the Earth limb. Gamma rays in the Earth limb are 
caused by cosmic-ray interaction in the Earth atmosphere and their spectrum is featureless, 
thus providing a useful control sample for spectral feature searches. While originally the 
limb signal was observed at similar signihcance as the feature in the Galactic center, the 
latest analysis of the Fermi-LAT yields somewhat lower signihcance. As the limb data 
can be observed for different event selection criteria, it can be used to calibrate the Monte 
Garlo-based estimate of the effective area. A feature might appear if the effective area 
is underestimated (or overestimated) in limited energy ranges, and indeed indications are 


reported that the efficiency is overestimated around the position of the line (Ackermann 


et al., 2013). However, the effect seems too small to explain the presence of the feature in 
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the limb data. 


As of the writing of this review, the origin of the line feature is not explained. However, its 
signihcance at the Galactic center is decreasing, as pointed out by Weniger[ (2013). Should 
an instrumental origin of the signal be ruled out (or results be inconclusive) independent 
conhrmation will be necessary. The currently operational lACT HESS H may be sensitive 


enough to soon provide such conhrmation (Bergstrom et al, 2012). The Fermi-LAT collab¬ 
oration has also presented a dedicated search for line features at energies between 100 MeV 


and 10 GeV (Albert et al, 2014). The analysis is particularly instructive as the number of 
expected events in that energy range is large and systematic uncertainties start to dominate. 


V. PERSPECTIVES 

A. Future instrumentation 

While many gamma ray features identihed by Fermi are intriguing, and in addition the 
recent measurements of the electron and positron spectra by Pamela, ATIG, Fermi and 
AMS-02 are sufficiently intriguing to generate a wealth of speculations, they also emphasize 
the need for improved energy resolution, position resolution, and background rejection up 
to the TeV range and even beyond. In this section, we discuss what future gamma-ray 
experiments will bring to the held of indirect dark matter detection, and the improvement 
that we can expect in our understanding of this held over the course of the next decade, and 
longer. 

1. Space Telescopes 

Over the next several years, the GALET and DAMPE experiments are expected to begin 
taking data. The GALorimetric Electron Telescope®, a Japan-led project that involves 
Italian and American institutes, is planned for launch in 2015 and will be installed on the 
Japanese Experiment Module on board the International Space Station. DAMPE (DArk 
Matter Particle Explorer, formerly known as TANSUO)®, is one of the hve satellite missions 

® http://calet.phys.Isu.edu/index.php 
® http://dpnc.unige.ch/dampe/index.html 
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selected by Chinese SPRPSS/CAS program, and is scheduled for launch slightly later than 
CALET, in 2015-2016. Italy and Switzerland take part in the DAMPE collaboration. Both 
projects feature a deep calorimeter to reach a total of 30 to 33 radiation lengths, in order 
to provide excellent energy resolution (better than 3% above 100 GeV) and electron/proton 
separation (~ 10^ rejection power) in the energy range of interest (1 GeV to 10 TeV). 
Both instruments will have roughly comparable performance, with a slightly larger electron 
geometrical factor for DAMPE (0.3 m^sr) than GALET (0.12 m^sr). Further information 


on the status of each project can be found in recent conference proceedings (Li and Yuan 


2012 Wu and Ghang, 2013; Yoshida, 2013) 


On a longer time scale (2018 and beyond), two other experiments will probe the high- 
energy electromagnetic sky: GAMMA-400 and HERD. Gamma-400^ is a Russian-led satel¬ 
lite observatory, planned for launch in 2018-2019. Building upon the successes of Fermi and 
AGILE, it features a gamma ray telescope reminiscent of the LAT, supplemented with a 
Konus-FG gamma ray burst monitor. The baseline design covers the range from 100 MeV 
to 10 TeV and is optimized for best performance around 100 GeV, where a very deep elec¬ 
tromagnetic calorimeter (25 radiation lengths compared to ~ 8.5 for the LAT), associated 
with a silicon strip tracker, will provide excellent energy and angular resolutions at such 
energies (a factor ten better angular resolution at 100 GeV than either DAMPE or GALET, 
and comparable energy resolution). Among the various gamma ray and cosmic ray science 
topics, dark matter searches, and especially the hunt for gamma ray lines, are a prime focus 


of the science case (Galper et al, 2013a Moiseev et al., 2013). Further information on the 


GAMMA-400 design and science case can be found in Galper et al. ( 2013a|[b ). 

The High Energy cosmic Radiation Detection (HERD) is an observatory planned for 
deployment on board the future Ghina space station ®. Design studies are still at an early 
stage of development, though the two primary science goals are already dehned as the search 
for a dark matter signal and the origin of Galactic cosmic rays. At this time, this mission 
is not in competition with any other Ghinese project. A recent development moved the 
baseline concept closer to the GAMMA-400 design, with a massive 3-dimensional calorimeter 
covered on hve sides by tracker silicon planes (more details can be found on the second HERD 


^ http: //gaiimia400 . lebedev. ru/indexeng. html 
® http://english.ihep.cas.cn/rs/fs/sm/SM/SM_aboutherd/ 
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international collaboration meeting, see http: //indico. ihep. ac. cn/conf erenceDisplay. 
py?confId=3808). 


Another experiment worth mentioning is PANGU (Wu et al., 2014). PANGU (the PAir- 
productioN Gamma-ray Unit) is a small mission optimized for spectro-imaging, timing and 
polarization studies in gamma rays in the still poorly explored energy band from 10 MeV to 
a few GeV. The present design is a pair conversion telescope with detector resolution about 
a of factor 2 better than previous instruments and a pointing resolution of a factor 3 to 5 
better than Fermi-LAT. 

While all the future space experiments mentioned above involve the high-energy regime, 
one should not forget that the phase space is still open for a dark matter candidate in the 


MeV domain (often coined sterile neutrino or LDM for Light Dark Matter, see e.g. Essig 
et al ., 2013[). This notoriously difficult energy regime is the focus of current intense research 


and development efforts geared toward the launch of a general-purpose astrophysical MeV 
space observatory by 2025 (see for instance the AstroMeV site http://astromev.in2p3.fr/). 

2. Ground Telescopes 

Two future instruments are expected to dominate the landscape of ground instruments: 
GTA® and HAWG^°. GTA (Gherenkov Telescope Array) is the next-generation lAGT (see 
e.g. Acharya et al. 2013; [Actis et al. , 2011). The energy range of this array is envisaged 
to be from a few tens of GeV to hundreds of TeV and the sensitivity is expected to be 
improved by one order of magnitude relative to current lAGTs. The angular and energy 
resolution are expected to lie between 0.1 (0.05) deg and 25 (10)% at low (high) energies, 
respectively. Technologically this energy range and sensitivity is achieved by combining a 
large number of single lAGTs of different size. Gurrently, one of the baseline designs foresees 
four H.E.S.S. II size telescopes (~ 23 meter diameter), about 30 medium size telescopes (~ 
12 meter diameter) and 30 to 70 small size telescopes (~ 7 meter). The US part of the 
consortium envisages to later extend the array with a large (~ 60) number of medium size 


telescopes, with particular view on high mass WIMPs (see e.g. Wood et al. (2013).) 


HAWG (High Altitude Water Gherenkov, Abeysekara et al., 2013a) is a second generation 


® https://portal.eta-observatory.org/Pages/Home.aspx 
http://www.hawc-observatory.org/ 
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cosmic ray and gamma ray observatory that builds upon the successful water Cherenkov 


technique pioneered by Milagro (Atkins et al, 2003). HAWC is designed to continuously 
(~ ‘24/7’ duty-cycle) survey during 10 years the 100 GeV to 100 TeV sky with a 1.8 sr 
instantaneous £eld-of-view telescope consisting of 300 water tanks, each instrumented with 
4 photomultiplier tubes. HAWC is a joint Mexico-U.S.A. project located at 4100 meters on 
the flanks of the Sierra Negra volcano near Puebla, Mexico. It started science operation in 
August 2013 with about 100 tanks, and the full array has very recently been completed. 
While the HAWC astrophysical science case is very strong thanks to its synoptic surveying 
nature, it may also prove competitive for indirect dark matter searches for WIMP masses 
larger than about 1 TeV. 


B. Until 2018 


Indirect detection of dark matter with gamma rays is entering a pivotal period. Upcoming 
instruments will reach sensitivities starting to probe into the most relevant parameter space 
at least for the most generic WIMP models. The Fermi Large Area Telescope will collect 
data until potentially beyond 2018 and with an updated event selection, known as “Pass 8” 


(Atwood et al, 2012 Atwood et al, 2013). The most anticipated results from the updated 


combined dSph analysis have been recently published (Ackermann et al, 2015), and conhrm 


that the dSph constraints are in mild tension with a dark matter interpretation of the GeV 
excess. Figure shows a simplistic forecast of the LAT combined dSph constraints with 
10 years of Pass 8 data and 3 more dSphs than currently known. Such an increase in 
the number of dSphs is motivated by the anticipation of new ultra-faint discoveries in the 
southern hemisphere with, for instance, surveys such as DES^^ and LSST^^. Indeed the DES 
collaboration has completed its first year of observations, and has -at the same time as 


other groups- recently announced the discovery of 8 potential new dSphs (Bechtol et al 
2015)(Koposov et al, 2015)^^, supporting an optimistic view of the hnal number of dSphs 


that will eventually be available to a gamma-ray analysis. As a matter of fact, the Fermi- 


LAT and DES collaborations released upper limits on these new candidates (Drlica-Wagner 


http://www.darkenergysurvey.org/ 
http://www.lsst.org/lsst/ 
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see also the recent discovery by PanSTARRS (Laevens et al. 2015) 
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et al., 2015), showing the promise for gamma-ray indirect searches of future optical surveys in 


the southern hemisphere. Despite no signifcant excess detected by the Fermz-LAT analysis 


with Pass 8 data selection, Geringer-Sameth et al. (2015) claim evidence (based on a 2-3 
a excess) using the less sensitive Pass 7 event selection. To our knowledge, the various 
future experiments reviewed in Section V.A.l have not shown sensitivity curves that could 
be overlaid on Figure 

For lACTs, 2012 showed the hrst light for H.E.S.S. II, and early results have been pub¬ 
lished recently. The addition of the 5th telescope will push the threshold down to about 50 
GeV, with corresponding impact on derived dark matter constraints. World-leading limits 
on WIMP annihilation can be expected from the Galactic center search for masses above 


about 800 GeV. Such a limit has been estimated by Gonrad (2014), and is reported in Fig¬ 


ure as well. H.E.S.S. II will also likely provide new insight into the issue of the possible 


line emission discussed in Section IV.G Following another path, VERITAS aims to accu¬ 


mulate a total of 1000 hours on Segue 1 (Smith et al., 2013) until 2018, a program that 


may rival the H.E.S.S. H constraints, as can also be seen in Figure Finally, the HAWG 
collaboration presented early results for such a search in the direction of Segue 1 with only 
30 water tanks and 82.8 days of observation (HAWG-30), and projection with the full array 
(HAWG-300) and one year of observation ( Abeysekara et al\ |2013b ). Using a J-factor of 
7.710^® GeV“^ cm“® sr they obtain limits, shown on Figurej^ that could become competitive 
above 1 TeV, when taking into account that HAWG does not need to allocate dedicated time 
“on source’,’ contrary to imaging Gherenkov telescopes (the latter do have a much better 
angular and energy resolution, though). Abeysekara et a/(] (2014) recently investigated more 
thoroughly HAWG sensitivity to a dark matter signal, based on simulations. 


C. Beyond 2018: future satellites and CTA 


Future satellite missions have been discussed in Section V.A.l GAMMA-400/DAMPE 


will be mostly contributing to searches for spectral features due to their superior energy 
resolution. HERD will provide a signihcant step forward with respect to Fermi -LAT also 
in terms of acceptance. For a dSph combined analysis, assuming the number of dSphs to 
be constant (which is a very conservative assumption), it seems that an exclusion of the 
thermal WIMP cross section should be possible up to WIMP masses of about 1 TeV by 
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2025. Detailed studies of the potential reach are however not known to us. This is not 
the case with CTA, for which the sensitivity reach has been studied in several publications, 


mostly focusing on the Galactic Center (Buckley et al., 2013 Doro et al., 2013 Pierre, 


Siegal-Gaskins, and Scott[ 2014 [Roszkowski, Sessolo, and Williams , 2015 Silverwood et al. 


2015 [Wood et al.\ [2013 ). dSphs and other approaches have been studied in Doro et al. 
( 2013[ ). Figure presents the limits anticipated at the Galactic center from Buckley et al. 
( 2013[ ). The one important caveat in interpreting the constraints is of course the systematics 
of the diffuse emission. As CTA lowers the threshold to about 10 GeV, constraints obtained 
from CTA will be prone to the same systematics as the Fermi-LAT Galactic center analysis. 


Silverwood et al. (2015) propose a quantitative discussion of this issue, and derive an upper 


limit curve that attempts to include such uncertainties, which is reproduced in Figure 



moM (GeV) 

FIG. 7 Selection of sensitivities to velocity averaged annihilation cross-section versus WIMP mass 
reachable with current or future experiments. The blue curve is from a Galactic center analysis 
with CTA, from Buckley et a/.j (2013). The green and black curves come from similar analyses, but 
attempting to account for the degradation of the limits due to uncertainties in the Galactic diffuse 


emission (Lefranc et al., 2015 Silverwood et al., 2015). The red curve corresponds to one-year of 


observation of Segue 1 with the full HAWC instrument (Abeysekara et al., 2013b). The cyan curve 


shows the Fermi-LAT Pass 8 limits by 2018, accounting for more dwarf spheroidals, based on a 
preliminary analysis of 5 years of data. The magenta curve shows the VERITAS expectations from 


1000 hours on Segue 1 (Smith et al, 2013), and the yellow curve is an estimate by Conrad (2014) 


of the limits H.E.S.S. II could reach on the Galactic center. 
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D. Progress on nuisance parameters 

As emphasized throughout this article, determination of particle dark matter limits and 
extraction of a signal in the future requires an understanding of astrophysical systematics. 
These systematics include measurements of dark matter distributions in the different targets 
and measurements of the diffuse gamma-ray backgrounds. In this subsection, we review the 
progress that is expected to be made with these systematics over the next several years. 


1. dSph mass distributions 


Much of this article has emphasized the methods of setting limits on the dark matter 
annihilation cross section using dSphs. These limits are constrained by the systematic 
uncertainties in the measurements of the J-factors, or the dark matter mass distributions in 
the dSphs. In this analysis above, we highlighted how this systematic uncertainty is treated 
as a statistical uncertainty in the gamma-ray analysis. 

In the future it is optimistic to expect that this systematic uncertainty in the J-factor will 
be reduced. For the brightest, classical dSphs, several thousands of line-of-sight velocities 
have now been measured. Order-of-magnitude improvements on these results will require 30- 
meter class telescopes which are expected to come online within the next decade. These 
larger samples of resolved stars with kinematic information will continue to improve the 
measurements of their dark matter distributions. A larger sample of line-of-sight velocities 
will be important to further test the initial results that have been obtained from multiple 
populations of stars in a small number of dSphs. 

A larger sample of line-of-sight velocities is important because forthcoming larger scale 
lACTs with angular resolution < 0.1° will be very sensitive to the dark matter distribution 
in the center of the dSphs. This is the regime in which the shape of the central dark 


matter prohle begins to signihcantly affect the determinations of J-factors (Charbonnier 


et al, 2011 Strigari, 2013). In addition to the improvement in the kinematics, it will also 


be important to improve the photometric data in the center of the dSphs. This is because 
there is a signihcant degeneracy between the central density of the stellar prohle and the 


central density of the dark matter prohle (Strigari, Frenk, and White, 2010). 
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http://www.gmto.org/, 

http: / / www.tmt.org/ 


http://www.eelt.org.uk/ 
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In addition to the increase in the sample of line-of-sight velocities from known dSphs, new 
ultra-faint satellites that will be discovered in future surveys will require kinematic follow 
up. It is intriguing to note that a forthcoming survey, such as the DES or LSST, may detect 
an object that is as massive and nearby as Segue 1. Beyond measurements of line-of-sight 
velocities, in the future it will also be possible to obtain transverse velocities of bright stars 
in several dSphs. This will go a long way towards breaking the degeneracy between the 
velocity anisotropy and the dark matter mass prohles. Adaptive optics systems on 30-meter 
class telescopes may be suited for these measurements. 


2. Local dark matter density distribution 


Gaia, an astrometric mission launched in 2013, will have an astrometric accuracy of 
approximately 1 mas and a photometric accuracy of 60 mmag for stars brighter than 20th 


magnitude (Perryman et al., 2001). Because of the improvement in the measurements of 


the phase space distribution of local stars, Gaia will by extension improve upon current 


measurements of the dark matter distribution in the Milky Way (Bovy and Rix, 2013), 


and especially the local dark matter density (Read, 2014). Measurements along these lines 


will be especially important for interpretation of the annihilation cross section limits from 
diffuse gamma-ray measurements, in addition to its obvious importance for local dark matter 
searches. 


3. Backgrounds 


As alluded to in Section III.A lAGTs and space borne instruments suffer from quite 
distinctive background. Until recently, the former could rely on a small held of view and 
steeply decreasing gamma-ray diffuse spectra to model the background below any source 
analysis as arising dominantly from an isotropic charged cosmic-ray component. With the 


detection of the Galactic diffuse emission by H.E.S.S. II (Abramowski et al., 2014c), in the 


near future this may be no longer true, at least for analyses at low Galactic latitude. 

On the contrary, cosmic-ray charged background is effectively rejected by an instrument 
like the LAT, and plays a potential role only at very high latitude, where the isotropic 
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component of a standard LAT sky model may dominate or for analyses related to the 


extragalactic gamma-ray background, especially at low energy (see Section IV.F.2). On the 
other hand, the gamma ray sky as seen by the LAT is extremely complex, and its imperfect 
angular resolution often demands a model over a much larger fraction of the sky than the 
immediate vicinity of a source of interest. 

While un-modeled sources above threshold can seriously alter the results of a dark matter 


search (Han et al, 2012) and unresolved sources can bias its signihcance, the most difficult 
component to account for remains the Galactic diffuse background. The uncertainties that 
enter its derivation are mentioned in Section IIII.DI and detailed in the review on indirect 


searches by Porter, Johnson, and Graham (2011). The effect on LAT all-sky residuals of 


variations in some of the components that define the model has been studied in [Ackermann 


et al. (2012a), while a Bayesian investigation of the posterior probability distribution of 


cosmic-ray-related parameters is the subject of Trotta et al. (2011). Whether predicted with 
ab initio propagation codes like GALPROP or estimated with template htting procedures 
as for the standard diffuse model of the LAT collaboration. Galactic diffuse models always 
leave residuals that need dedicated treatment. Nevertheless, these residuals are sufficiently 
low on average to guarantee that a nominal diffuse model is useful for the large majority of 
LAT analyses. What remains very difficult to achieve is proper propagation of uncertainties. 


which often result from an ad hoc study, as in the EGB determination (Ackermann et al. 


2015a), the construction of a catalog of supernova remnants (Hewitt et al., 2013), or the 


extraction of upper limits on a dark matter signal from the EGB (Ackermann et al ., 2015c 


Ajello et al ., 2015 Di Mauro and Donato, 2015). In the future, new insight may come from 


next-generation propagation codes, like PIGARD ( Kissmaimf 2014), or new cosmic-ray data, 
notably from AMS-02^®. Also improved constraints are expected on the interstellar radiation 
held, a critical ingredient to the determination of the large-scale inverse Gompton emission 
of the Galaxy, coupled to further investigations on statistical ways to properly account for 
these uncertainties. 


this isotropic component is thus derived by the LAT collaboration for major event classes, and distributed 
at the FSSC : http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
http://www.ams02.org/ 
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E. Beyond 2024 


Beyond the next generation satellite experiments (DAMPE, HERD, GAMMA-400, 
PANGU etc.) and GTA, we are not aware of concrete proposals for next to next gener¬ 
ation gamma ray telescopes for dark matter detection. Indeed, the question can be asked 
if a “Super-GTA” or “Dark Matter Array,” is possible. This was for example proposed 
in Bergstrom, Bringmann, and Edsjo (2011), where GTA’s sensitivity was improved by a 


factor 10 as a “Gedanken-experiment.” In many ways, GTA will be the ultimate indirect 
detection experiment. Progress signihcantly below the thermal cross-section is hampered 
by the irreducible electron background and the systematic uncertainties in the knowledge of 
the acceptance. 

VI. CONCLUSIONS 


In this article we have reviewed the methods that are utilized to search for gamma 
rays from dark matter annihilation, discussed the challenges in extracting the signal, and 
reviewed the results that have been obtained by different experiments. We have in particular 
focused on the great deal of progress over the past few years that has resulted from the 
analyses of FermALAT data. Without question, during the course of its over 6 year mission 
lifetime, the FermALAT has been a phenomenal success, breaking new ground in indirect 


dark matter searches and far exceeding the pre-launch projections for its reach (Baltz et al. 


2008). Analysis of several systems, most notably dwarf spheroidals, are now providing robust 


upper limits on WIMPs in the mass range ~ 10 — 100 GeV, ruling out velocity independent 
annihilation cross sections near the cosmologically-motivated thermal relic cross section scale 
of {av) ~ 3 X 10“^® cm^ s“h 

Analysis of FermALAT data has not, however, been without controversy. After over 
6 and a half years of science operations, probably the biggest debate revolves around the 
analysis of data from the inner Galaxy, in particular the nature of the possible extended 
emission observed at a few GeV. When interpreting this emission as due to dark matter, 
at present it is possible to £t this extended excess to a wide range of WIMP masses, cross 


sections, and annihilation channels (Galore et al, 2014). Given this wide range of models 


that seemingly work to £t this emission, as well as the uncertainties in Milky Way dark 
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matter distribution and astrophysical sources of gamma rays in the Galactic center, it is 
clear that establishing this result as due to dark matter will require a conhrmation from 
other sources. The dwarf spheroidals will provide an especially important conhrmation, and 
indeed the current limits from dwarf spheroidals are able to rule out regions of parameter 
space that explain the extended emission. 

Even after over six years of tremendously successful data taking, there is still much 
to be gained from future Fermz-LAT data. WlMPs with mass in the range 10-100 GeV 
with velocity-independent annihilation cross sections are now strongly ruled out by dwarf 
spheroidal analysis. WlMPs with mass greater than about 100 GeV have been more difficult 
to constrain because of the photon counting statistics above about a few GeV. In this regime 
the improvement in sensitivity will scale linearly with time. Forthcoming FermGLAT data 
will thus prove vital in extending WIMP limits to this higher mass regime. lAGTs main 
target is the Galactic center and analyses here might be systematics limited at masses below 
1 TeV, and thus LAT data with a 10 year exposure will be competitive up to a WIMP mass 
of about 1 TeV. 

Over the entire detected energy regime, we have repeatedly seen that potential dark 
matter signals in the FermGLAT data have been compromised by unresolved point sources 
that are near the threshold limit for point source detection. Because of the foreground 
model, there is a flux threshold below which it will not be possible to identify individual 
point sources, as point sources below this threshold contribute to the diffuse foreground. 
Future Fermz-LAT data will be important for better identifying faint point sources, and 
thereby understanding their contribution to the diffuse Galactic model. 

Turning to lAGTs, as already mentioned the mass range above a few TeV will be covered 
robustly and with unrivaled sensitivity by GTA. This will be especially important if the 
LHG searches do not hnd hints of new physics, thereby pushing potential WIMP masses to 
above 500 GeV to 1 TeV. 

Taking an even bigger step back, in order to verify any signal that is obtained from 
indirect detection experiments, it will be crucial to examine the signal in the context of 
other dark matter searches, in particular those from direct detection and colliders. While 


the detailed discussion of results of direct searches (Gushman et al, 2013) and accelerator 


searches (Dutta, 2014) -as well as other indirect probes- is beyond the scope of this article, 
it is important to ultimately emphasize the complementarity of the different methods. 
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Though a comparison between the sensitivities of the different methods is inherently 
model dependent, we can highlight it in two different theoretical set-ups. For example 
in minimal supersymmetry, there is little correlation between scattering cross section and 


annihilation cross-section (Bergstrom, 2012). Direct detection searches are seen to constrain 
the model space approximately orthogonal to indirect detection by gamma rays. Extending 
beyond direct and indirect searches, collider detection and subsequent measurements of the 
sparticle mass spectrum and splittings of a supersymmetric model providing dark matter 


might be used to calculate annihilation cross-sections and relic density (Baltz et al ., 2006). In 
the framework of phenomenological minimal supersymmetry, however, it is quite conceivable 
that only after a combination with indirect detection experiments the solution providing dark 
matter can be identified (Bertone et al.\ 2012). 
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